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of  the  data  acquisition  system;  and  a  modified,  more  rapid  data  reduction  program. 
Admittances  and  bulk  loss  coefficients  of  nonaiuminized  and  aluminized  propellants 
were  measured  over  a  range  of  frequencies.^ 

<4TASK  II.  Studies  were  continued  on  external  burning  using  an  axisymmetric  model 
at  VAacb  3.  i  Early  work  using  contoured  test  section  walls  and  cold  gas  base  region 
injection  is  reviewe(Jtord6rff5igtrate  wake  detail  and  lengthscalc  changes  unarTT) 
the  influence  of  simulated  external/base  burning j^fTests  with  actual  combustion  of 
radially  and  axially  injected  hydrogen  are  then  reported.  Outstanding  performance 
values  with  significant  base  drag  reduction  is  shown  for  injection  and  burning  directly 
in  the  near-wake  (base  burning). .  Current  at  temp  ts  at  radial  injection  and  burning 
in  the  free  stream  (external  burningDhave  not  yet  succeeded.  These  tests,  however, 
have  defined  an  envelope  within  whim  external  burning,  if  feasible,  will  presently 
be  achieved.  _ _ ^ 

TASK  III.  Investigations  continued  of  the  accumulation  processes  of  aluminum  on  the 
propellant  burning  surface  that  lead  to  formation  of  agglomerate  droplets,  the  size 
of  which  dominates  aluminum  combustion..  Using  dry-pressed  mixtures  of  ammonium 
/^perchlorate,  aluminum  and  carnauba  wax  powders,  the  effects  of  particle  sizes, 

/  mixture  ratios  and  pressure  were  studied  by  hign  speed  photography.  Results  were 
I  correlated  with  a  scenario  describing  the  effects  of  propellant  microstructure, 

I  retention  processes  for  aluminum  on  the  burning  surface,  response  of  aluminum  powders 
I  to  heating,  and  inf  tarnation  process  for  sintered  aluminum. 

f  Quench  tests  at  atmospheric  pressure  on  the  combustion  plume  above  the  burning 
[  surface  showed  that  the  agglomerates  in  these  tests  burned  much  as  aluminum 
I  droplets  in  earlier  single  particle-controlled  atmosphere  tests.  An  exception  resulted 

1  in  tests  with  large  (95  urn)  aluminum  particles  which  did  not  ignite  completely.  An 

\  apparatus  for  tests  at  elevated  pressure  was  designed  and  partially  fabricated. 

^  TASK  IV.  Experimental  and  analytical  studies'  were  continued  on  the  subject  of 
turbulence-induced  pressure  fluctuations  in  a  roc^-f*1"*  ‘Tn  overcome  prior 

deficiencies'  in  theory  ^£,new  theory  "Was  constructed  for  the  generation  of  pressure 
disturbances  by  the  turbUfence.  Undeniable  separation  has  been  achieved  between 
propagational  and  local  pressure  fluctuations,  and  the  agreement  between  theory 
and  experiment  are  adequate.  Measurements  include  microphone  and  hot  film 
anemometry,  and  time  serietfbnalysis  is  applied  to  the  results. 
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General  Introduction 


Progress  is  reported  for  the  second  year  effort  of  four  distinct  projects  under 
contract  AFOSR  F49620-78-C-0003.  The  individual  projects  are  identified  as 
distinct  tasks|.this  distinction  being  preserved  because  it  reflects  also  distinct  test 
facilities  and  personnel  assignments.  Each  task  report  contains  its  own  abstract, 
introduction,  and  progress  report. 

Task  1  is  concerned  with  combustion  stability,  with  particular  reference  to 
measurement  of  the  response  of  solid  propellant  combustion  to  oscillatory  flow 
environments.  Primary  effort  is  on  development  of  the  impedance  tube  method  for 
measuring  response  functions.  This  includes  improving  its  accuracy,  economy  of 
operation,  and  versatility  of  application.  Exploration  of  other  methods  of 
measuring  response  functions  is  also  conducted. 

Task  n  is  concerned  with  combined  external  and  base  burning  as  a  means  of 
developing  sustained  thrust  for  missiles  with  minimum  weight  and  volume  demands 
on  the  propulsion  system.  A  facility  was  developed  previously  for  injection  and 
combustion  of  hydrogen  in  and  around  the  near-wake  shear  flow  of  an  axisymmetric 
body.  Initial  studies  were  made  of  the  flow  field  using  nonreactive  gases,  and  the 
present  report  covers  recent  studies  with  hydrogen  injection  and  burning. 

Task  HI  is  concerned  with  the  roie  of  the  aluminum  powder  ingredient  in 
propellant  combustion  and  2-phase  product  flow.  The  steps  in  the  metal 
"metabolism"  are  traced  through  the  phase  of  accumulation  on  the  burning  surface, 
the  melting  and  coalescing  of  particles  into  relatively  large  droplets,  the  ignition 
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and  combustion  of  the  droplets,  effect  of  combustor  gas  flow  on  these  steps,  and 
their  response  to  gas  oscillations.  The  present  report  summarizes  combustion 
photographs,  tests  of  burning  surface  behavior  of  model  propellants  (dry-pressed 
powders)  involving  various  combinations  of  aluminum,  ammonium  perchlorate  and 
hydrocarbon  wax  powder.  Also  reported  are  results  of  quench  tests,  including 
quench-coilection  of  condensed  material  in  the  combustion  phase. 

Task  IV  is  concerned  wtih  the  source  of  pressure  fluctuations  in  rocket 
motors.  It  is  believed  that  a  significant  contribution  is  made  by  the  interaction  of 
turbulence  of  the  flow  with  the  nozzle.  A  facility  was  constructed  previously  to 
measure  this  effect,  and  modifications  are  continuing.  Recent  testing  is  reported 
here,  and  a  new  theory  is  presented. 
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Abstract 

Results  obtained  under  an  AFOSR  sponsored  research  program  that  is  concerned 
with  the  development  of  experimental  techniques  for  the  determination  of  the 
admittances  of  burning  solid  propellants  under  conditions  simulating  those  observed  in 
unstable  rocket  motors  and  the  determination  of  the  admittances  of  various  solid 
propellant  are  described.  As  part  of  this  program,  much  effort  has  gone  into  improving 
the  accuracy  of  the  impedance  tube  technique  that  was  developed  earlier  under  this 
program.  These  improvements  included  increasing  the  number  of  transducers  utilized  to 
measure  the  standing  wave  structure  to  fifteen,  increasing  the  dynamic  measurement 
range  to  cover  110-165  db,  increasing  the  sampling  rate  of  the  mini  .computer  based  data 
acquisition  system,  and  reducing  the  number  of  numerical  iterations  required  for 
determining  the  propellant  admittance  from  the  measured  acoustic  pressure  data.  The 
improved  impedance  tube  set-up  was  then  used  to  determine  the  frequency  dependence  of 
the  admittances  and  bulk  loss  coefficients  of  aluminized  and  non-aluminized  propellants 
and  some  of  these  results  are  presented  in  this  report.  Finally,  efforts  dealing  with  the 
development  of  a  direct  admittance  measurement  technique  are  described. 


1.  Introduction 


This  report  summarizes  progress  made  under  Task  I  of  AFOSR  Contract  No. 
F49620-78-C-0003  during  the  period  October  1,  1978  to  September  30,  1979.  This  research 
was  concerned  with  die  development  of  the  experimental  techniques  for  the  determination 
of  the  admittances  of  burning  solid  propellants  under  a  variety  of  conditions  simulating 
those  observed  in  unstable  solid  rockets  and  the  measurement  of  the  admittances  of 
various  aluminized  and  non-aluminized  composite  solid  propellants.  The  ability  to  predict 
and  measure  the  burning  rate  characteristics  of  solid  propellants  under  a  variety  of 
unsteady  rocket  motor  operating  conditions  is  of  great  scientific  and  practical 
importance.  Such  data  is  needed  for  the  apriori  prediction  of  the  combustion  stability 
characteristics  of  proposed  solid  propellant  rocket  motor  designs  as  well  as  the 
determination  of  the  relative  driving  capabilities  of  various  solid  propellant  formulations. 
Also,  there  is  a  need  to  correlate  the  predictions  of  theoretical  solid  propellant 
combustion  models  and  relevant  experimental  data.  Establishing  the  validity  of  theoretical 
combustion  models  would  demonstrate  a  level  of  technology  which  would  greatly  redude 
future  development  testing  requirements.  The  progress  made  in  reaching  the  goals  of  this 
program  are  briefly  outlined  in  the  remainder  of  this  report. 

As  the  objectives  and  background  for  this  research  program  have  been  discussed  in 
much  detail  in  Ref.  1  only  a  brief  review  of  this  material  will  be  provided  herein.  The 
propellant  in  unstable  solid  rockets  is  exposed  to  regions  of  primarily  pressure  oscillations, 
regions  of  primarily  velocity  oscillations  and  regions  that  contain  both  pressure  and 
velocity  oscillations.  To  predict  solid  rocket  motor  stability,  data  describing  the  propellants 
burning  rate  at  each  of  these  regions,  at  different  frequencies,  different  oscillation 
amplitudes  and  different  engine  operating  conditions  are  needed.  Most  of  the  research 
efforts  in  this  area  to  date  have  concentrated  on  the  pressure  coupled  response  that 
describes  the  propellant  burning  rate  response  when  exposed  to  small-amplitude  pressure 
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oscillations  only.  Most  of  the  experimental  efforts  in  this  area  to  date  depend  upon  the 

2 

utilization  of  the  T-burner  which  provided  much  useful  data  in  spite  of  some  obvious 
shortcomings*.  Careful  considerations  of  what  has  been  done  in  this  area  to  date  clearly 
indicate  the  need  for  the  development  of  alternate  experimental  approaches  that  are 
hopefully  devoid  of  some  of  the  difficulties  that  affect  T-burners.  Results  provided  by 
these  new  measurement  techniques  could  also  be  used  to  check  the  accuracy  of  the  T- 
burner  data. 

The  work  conducted  under  this  program  can  be  generally  divided  into  two  parts; 
namely,  (1)  the  development  of  experimental  techniques  for  the  measurement  of  the 
admittances  of  solid  propellants  under  different  conditions  simulating  those  observed  in 
unstable  rockets  and  (2)  the  acquisition  of  solid  propellant  admittance  data. 

The  efforts  concerned  with  the  development  of  new  experimental  techniques 
included  the  improvement  of  the  impedance  tube  technique  developed  earlier  under  this 
program  and  the  investigation  of  the  so-called  "direct"  burn  rate  measurement  techniques. 
In  this  case  the  term  "direct"  implies  experimental  techniques  that  are  capable  of  directly 
measuring  the  mass  loss  of  the  propellant  sample  or  its  surface  regression.  These  efforts 
are  still  in  progress  and  their  state  of  development  is  described  in  the  following  section  of 
this  report. 

The  efforts  to  acquire  solid  propellant  admittance  data  concentrated  on  the 
determination  of  the  admittances  of  two  non-aJuminized  and  two  aluminized  propellants 
at  300  psig  over  the  frequency  range  from  300  to  1000  Hz.  Of  special  significance  are  the 
investigations  of  the  admittances  of  aluminized  propellants  where  the  utilization  of  the 
impedance  tube  allowed  the  simultaneous  determination  of  the  gas  phase  losses  caused  by 
the  aluminum  oxide  particles  in  the  impedance  tube  and  the  admittance  at  the  propellant 
surface.  The  capability  to  perform  these  two  measurements  simultaneously  offers  the 
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possibility  for  the  determination  of  the  manner  in  which  the  addition  of  aluminum  or 
metai  particles  affects  combustion  instabilities;  that  is,  whether  the  addition  of  aluminum 
particles  attenuate  or  amplify  combustion  instabilities  by  providing  gas  phase  damping  or 
amplification  and/or  modifying  the  combustion  response  of  the  solid  propellants.  The 
resolution  of  this  problem  will  require  the  utilization  of  the  developed  impedance  tube 
technique  in  a  series  of  tests  in  which  the  same  propellants  will  be  tested  with  and  without 
metal  additives  and  the  results  will  be  compared.  Some  of  the  propellant  admittance  data 
measured  in  the  course  of  this  investigation  will  be  presented  in  the  following  section. 

The  ability  to  utilize  the  impedance  tube  technique  in  the  measurement  of  the 
admittances  of  aluminized  propellants  is  of  special  significance  as  other  experimental 
techniques,  such  as  the  pulse  and  variable  area  T-burners,  have  experienced  considerable 
difficulties  in  the  determination  of  the  admittances  of  aluminum  solid  propellants. 

Results  obtained  under  this  program  were  presented  during  the  reporting  period  at 
the  1 7th  Aerospace  Sciences  Meeting  in  New  Orleans,  La.  and  the  16th  JANNAF 
Combustion  Meeting  at  Naval  Post  Graduate  School  in  Monterey,  Calif.  These  works  also 
appeared  as  the  following  publications: 

1.  M.  Salikuddin  and  B.  T.  Zinn," Adaptation  of  the  Impedance  Tube  Technique  for  the 
Measuement  of  Combustion  Process  Admittances,"  A1AA  Paper  No.  79-0167, 
January  1978. 

2.  "  Experimental  Determination  of  Solid  Propellant  Admittances  by  the  Impedance 
Tube  Method',  to  appear  in  the  Proceedings  of  the  16th  JANNAF  Combustion 
Conference,  Sept.  10-14,  1979. 

Paper  No.  1  above  has  also  been  accepted  for  publication  in  the  December  issue  of 
the  Journal  of  Sound  and  Vibrations. 
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In  addition,  the  paper 

3.  Baum,  3.  D.,  Daniel,  B.  R.  and  Zinn,  B.  T„  "Determination  of  Solid  Propellant 
Admittances  of  Impedance  Tube  Method,"  has  been  accepted  for  presentation  at 
thelSth  AIAA  Aerospace  Sciences  Meeting  that  will  be  held  in  Pasadena, 
California. 

All  of  the  above  papers  are  included  in  Appendix  1-A  of  this  report. 

Finally,  the  Ph.D.  thesis  of  Dr.  M.  Salikuddin  that  was  supported  under  this  program 
was  published  and  approved  by  the  Georgia  Tech  faculty  in  April  1978.  This  thesis  is 
entitled: 

"Application  of  the  Impedance  Tube  Technique  in  the  Measurement  of  Burning  Solid 
Propellant  Admittances?' 

In  the  following  section  a  brief  summary  of  the  progress  made  during  the  reporting 
period  is  provided.  More  detailed  descriptions  of  these  efforts  can  be  found  in  the  papers 
included  in  Appendix  I-A. 


(A)  Improvements  in  the  Impedance  Tube  Technique 

The  principles  of  the  operation  of  the  impedance  tube  technique  are  provided  in 
Paper  No.  1  in  App.  I-A.  As  with  any  newly  developed  experimental  technique,  much 
effort  was  expended  during  the  reporting  period  on  the  improvement  of  the  operation  and 
accuracy  of  the  impedance  tube  technique.  These  efforts  can  be  divided  into  three  parts; 
namely,  (a)  improvement  of  the  experimental  procedure;  (b)  improvement  of  the  data 
acquistion  procedure;  and  (c)  improvement  of  the  data  reduction  procedure.  Item  (a) 
involved  the  use  of  additional  transducers  (up  to  15)  in  the  acoustic  pressure 
measurements  and  the  development  of  a  new  calibration  procedure  that  increased  the 
previously  available  dynamic  measurement  range  to  cover  the  range  of  110-165  db.  Item 
(b)  involved  the  inclusion  of  a  "direct  memory  access"  into  the  mini  computer  based  data 
acquisition  system  which  greatly  increased  the  data  sampling  rate.  This  change  resulted  in 
considerable  improvement  in  the  quality  of  the  measured  acoustic  data.  Finally,  Item  (c) 
involved  the  mathematical  reformulation  of  the  data  reduction  procedure  in  order  to 
eliminate  the  majority  of  the  previously  used  numerical  iterations.  This  improvement 
considerably  reduces  the  amount  of  computer  time  required  to  obtain  the  desired 
admittances  from  the  measured  acoustic  pressure  data. 

Since  the  efforts  covered  under  items  (a),  (b)  and  (c)  above  are  described  in  detail 
in  Papers  2  and  3  of  App.  I-A  of  this  write  up,  further  discussion  of  these  efforts  will 
not  be  repeated  herein.  In  addition,  these  references  contain  data  illustrating  the 
repeatability  and  consistency  of  the  data  measured  by  the  impedance  tube  technique. 

(B)  Measured  Solid  Propellant  Admittance  Data 

During  the  report  period,  the  impedance  tube  technique  was  utilized  to  measure  the 
admittances  of  two  non-aluminized  composite  propellants  (i.e.,  A- 13  and  A- 14)  provided  to 
this  group  by  NWC,  China  Lake,  Calif,  and  two  aluminized  propeiiants.  The  latter 
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consisted  of  a  UTP-3001  propellant  that  is  a  Titan  III  C  Booster  Propellant  provided  by 
Prof.  E.  W.  Price  at  Georgia  Tech  who  uses  it  in  his  research  program,  and  a  UTP- 19360 
propellant  provided  by  Dr.  L.  Strand  from  the  3PL  Laboratories  in  Pasadena,  California. 

Typical  results  obtained  with  the  UTP-  3001  aluminized  propellant  are  presented  in 
Figs.  1  through  4  of  this  report.  In  this  case,  Figs.  I  and  2  present  a  comparison  between 
the  instantaneous  pressure  amplitude  and  phase  profiles  predicted  utilizing  the  determined 
propellant  admittances  and  the  measued  experimental  data.  Excellent  agreement  between 
the  two  sets  of  data  is  demonstrated.  Furthermore,  attention  should  be  given  to  the 
increase  in  the  magnitude  of  the  second  amplitude  minima,  as  compared  to  the  first 
amplitude  minima.  This  change  was  caused  by  the  gas  phase  attenuation  provided  by  the 
presence  of  aluminum  oxide  particles  in  the  flow.  The  large  gas  phase  damping  is  also 
supported  by  the  change  in  the  positive  slope  of  the  phase  space  curve  near  the  respective 
minimum  points.  The  frequency  dependence  of  the  UTP-3001  propellant  driving 
characteristics,  as  described  by  the  real  part  of  the  admittance  YR,  and  the  gas  phase 
damping,  as  described  by  the  gas  phase  damping  coefficient  G,  are  plotted  in  Figs.  3  and  4, 
respectively.  Similar  data  for  the  UTP- 19 360  propellant  has  also  been  obtained  during  the 
research  period  covered  by  this  report. 

(C)  Evaluation  of  the  Feasibility  of  the  "Direct”  Impedance  Measurement  Techniques. 

As  stated  earlier,  the  "direct"  measurement  techniques  involve  the  direct 
measurement  of  the  propellant  mass  loss  and/or  the  regression  of  its  surface,  when 
exposed  to  an  oscillatory  flow  field.  The  difficulty  in  performing  these  measurements  is 
that  they  require  the  measurement  of  extremely  small  mass  or  distance  changes  over 
extremely  short  time  periods.  In  contrast  to  the  "direct"  measurements,  most  of  the 
currently  used  techniques  for  admittance  measurement  (e.g.,  the  T-burner  and  the 
impedance  tube)  utilize  measured  acoustic  pressure  data  to  determine  the  propellant 
admittance.  In  these  cases  the  accuracy  of  the  determined  admittances  depends,  among 
other  things,  upon  the  accuracy  of  the  theoretical  model  used  to  describe  the  system 
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acoustics. 

Some  of  tiie  research  efforts  conducted  under  this  program  were  concerned  with 
the  investigation  of  the  feasibility  of  the  some  potentially  promising  direct  measurement 
techniques.  In  the  first  part  of  this  program  the  feasibility  of  utilizing  the  Optron,  an 
optical  target  follower,  to  follow  the  regression  of  a  burning  solid  propellant  surface  was 
considered.  The  characteristics  of  the  Optron  System  were  discussed  in  last  year's 
proposal3  and  they  will  not  be  repeated  herein. 

Part  of  the  efforts  conducted  under  this  part  of  the  program  involved  the  borrowing 
of  an  Optron  system  that  was  available  at  AFRPL.  Upon  receipt  of  this  system  much 
effort  went  into  its  checkout  and,  unfortunately,  this  system  did  not  perform  adequately. 
Simultaneously  with  the  Optron's  checkout,  the  ability  of  the  Optron  system  to  accurately 
measure  the  small  changes  in  propellant  surface  regression  in  response  to  an  oscillatory 
flow  field  were  considered.  This  analysis  indicated  that  the  distance  changes  that  one 
would  need  to  measure  in  order  to  determine  the  propellant  admittance  are  of  the  same 
order  of  magnitude  as  the  surface  irregularities  caused  by  the  presence  of  ammonium 
pechlorate  in  the  binder.  Consequently,  the  ability  of  the  Optron  to  distinguish  the 
"oscillatory"  changes  in  the  propellant  regression  caused  by  the  propellant  response  to  an 
oscillatory  environment  from  "random"  regression  rate  fluctuations  is  uncertain.  To 
answer  this  question,  a  newer  and  more  accurate  Optron  system  will  have  to  be  utilized. 
Considering  the  uncertainty  of  the  results  of  such  an  effort  and  the  high  cost  of  a  new 
Optron  system,  a  decision  has  been  made  to  terminate  these  efforts  for  the  time  being 
and,  instead,  focus  our  efforts  on  the  possibility  of  directly  measuring  propellants  mass 
losses. 

The  other  portion  of  this  study  was  concerned  with  the  development  of  an  accurate 

experimental  technique  for  direct  mass  loss  measurement.  Efforts  conducted  in  this  area 
4  5 

in  the  Soviet  Union  ’  concentrated  on  the  development  of  a  system  whose  mass  change 


could  be  related  to  a  frequency  change  that  can  be  measured  with  much  accuracy.  A 


no 


schematic  of  one  of  these  setups  where  the  change  in  the  natural  frequency  of  an 
equivalent  mass-spring-dashpot  system  is  measured  as  shown  in  Fig.  5.  In  this  experiment, 
a  change  in  the  mass  of  the  burning  propellant  is  reflected  by  a  proportional  change  in  the 
natural  frequency  of  the  system.  The  major  shortcoming  of  this  effort  is  the  low  natural 
frequency  of  the  developed  mechanical  system.  Therefore,  the  relatively  small  change  in 
the  natural  frequency  due  to  the  propellant  mass  flux  precludes  the  high  degree  of 
resolution  required  for  meaningful  propellant  admittance  measurements. 

The  most  promising  technique  for  direct  mass  loss  measurement  is  a  modification 
of  the  mass-spring-dashpot  system  which  uses  the  high  natural  frequency  characteristics 
of  piezoelectric  crystals.  A  schematic  of  the  setup  currently  being  investigated  under  this 
program  is  also  shown  in  Fig.  5  .  In  initial  investigations  with  this  setup,  quartz  crystals 
were  used  as  the  basic  oscillatory  system.  Unfortunately,  the  oscillatory  characteristics  of 
the  quartz  crystals  is  such  that  the  crystals  will  not  oscillate  with  the  relatively  high 
propellant  mass  loadings  that  are  necessary  for  mass  flux  measurements. 

More  recent  work  in  this  area  has  been  directed  toward  using  a  ceramic  crystal  as 
the  oscillating  system.  Progress  in  this  area  to-date  includes  the  selection  and  purchase  of 
several  configurations  of  piezo-ceramic  crystals  with  unloaded  natural  frequencies  of  the 
order  of  9  Mega  Hertz.  It  is  anticipated  that  at  this  relatively  high  frequency  the  change 
in  the  natural  frequency  of  the  system  due  to  propellant  mass  loading  will  be  of  an  order 
which  will  make  measurement  of  the  mass  flux  of  a  burning  solid  propellant  possible. 
However,  the  feasibility  of  the  application  of  the  ceramic  crystals  for  direct  mass 
measurement  will  depend  on  the  static  loading  tests  which  will  be  conducted  as  soon  as 
the  crystal  elements  are  delivered  from  the  manufacturer.  This  work  is  currently  in 


progress. 
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HI,  Future  Research 

This  program  is  still  committed  to  the  objective  of  improving  current  experimental 
capabilities  for  measuring  the  admittances  of  various  solid  propellants  under  a  variety  of 
conditions  simulating  these  observed  in  unstable  rocket  motors  and  the  acquisition  of 
propellant  response  factor  data  for  a  variety  of  aluminized  and  non-aluminized  propellants 
over  frequency  ranges  that  are  of  interest  to  the  propulsion  community.  Specifically,  the 
following  studies  will  continue  under  this  program:  (a)  the  investigation  of  the  feasibility 
of  utilizing  a  piezo-ceramic  crystal  for  the  direct  measurement  of  the  propellant  mass 
loss;  (b)  the  determination  of  the  admittances  of  various  aluminized  and  non-aluminized 
propellants  with  special  emphasis  given  to  the  determination  of  the  effect  of  metal 
addition  upon  the  propellant  admittance;  and  (c)  an  effort  will  be  undertaken  to 
determine  whether  the  impedance  tube  technique  can  be  used  to  determine  the  velocity- 
coupled  response  of  solid  propellants. 
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Figure  4.  The  Frequency  Dependence  of  Che  Bulk  Loss  Coefficient  G  for  a 
UTP-3001  Aluminized  Propellant  Tested  at  300  psig. 
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Abstract 


Experimental  studies  of  base  pressure  manipulation  for  an  axisymmetric 
model  at  Mach  3  with  simulated  and  actual  external  base  burning  are  described. 
Early  work  using  contoured  test  section  walls  and  cold  gas  base  region  injection  is 
reviewed  to  demonstrate  wake  detail  and  length  scale  changes  under  the  influence 
of  simulated  external/base  burning.  Tests  with  actual  combustion  of  radially  and 
axially  injected  hydrogen  are  then  reported.  Outstanding  performance  values  with 
significant  base  drag  reduction  is  shown  for  injection  and  burning  directly  in  the 
near-wake  (base  burning).  Current  attempts  at  radial  injection  and  burning  in  the 
free  stream  (external  burning)  have  not  yet  succeeded.  These  tests,  however,  have 
defined  an  envelope  within  which  external  burning,  if  feasible,  will  presently  be 
achieved. 


I.  Introduction 


The  application  of  burning,  either  in  the  near  wake  (base)  or  in  the  free 
stream  adjacent  to  the  near  wake  (external),  can  provide  significant  base  drag 
reduction  for  projectiles  operating  at  airbreathing  altitudes.  The  coupling  of 
attractive  performance  values  with  hardware  simplicity  makes  this  a  particularly 
interesting  propulsion  concept. 

The  feasibility  of  external/base  burning  for  propulsion  in  which  the  subsonic 
near-wake  transmits  elevated  downstream  pressures  to  the  projectile  base  has  been 
established.  At  present,  however,  no  analytical  model  exists  which  adequately 
treats  the  complex  base  flows  that  result  from  combustion  in  and  near  the  wake.  A 
major  contribution  in  understanding  these  flows,  establishing  and  optimizing 
performance  criteria  and  providing  direction  for  theory  can  be  supplied 
experimentally  through  detailed  pressure  and  temperature  measurements  in  the 
wake  region  for  various  modes  of  external/base  combustion.  Current  studies  in  the 
School  of  Aerospace  Engineering  at  Georgia  Tech  seek  to  provide  this  information. 

A  versatile  experimental  facility  has  been  established  for  both  simulation  and 
actual  combustion  testing  of  axisymmetric  external/base  burning  configurations  at 
Mach  3  with  a  5.72  cm  dia.  model  in  a  15.2  cm  dia.  test  section  (Reynolds  number 
based  on  model  diameter  =  3.0  x  10^).  The  blow-down  facility  is  shown  in  Figures  1 
and  2  in  Appendix  II-A.  The  hollow  cylindrical  model  is  mounted  in  the  subsonic 
flow,  passes  through  the  nozzle  throat  and  terminates  in  the  plexiglass  test  section. 
Gases  for  base  region  injection  are  brought  in  through  four  of  the  eight  support 
strut  access  ports.  The  remaining  four  ports  are  used  to  bring  model  surface  static 
pressure  tubes  to  the  external  data  acquisition  system.  Detailed  pressure  and 
temperature  wake  measurements  are  made  with  computer  controlled  traverses  of 
probes  introduced  through  the  constant-area  diffuser  downstream  of  the  test 


section.  Run  times  range  from  4-7  minutes  depending  on  the  stagnation  pressure 
selected.  Tunnel  flow  is  not  heated  and  stagnation  temperatures  drift  downward 
from  ambient  to  approximately  -20°C  during  typical  tests. 


11-5 


II.  Program  Progress 

Investigations  of  external/base  burning  with  hydrogen  have  continued  this 
year  in  the  Georgia  Tech  Mach  3  Test  Facility.  The  experiments  were  designed  to 
meet  the  following  goals: 

1.  Measurement  of  over-all  performance  and  wake  detail  with  pure  base 
burning 

2.  Determination  of  pure  external  burning  feasibility  with  subsequent 
performance  and  flow  detail  measurements 

3.  Investigation  of  base  burning  in  combination  with  external  burning  to 
define  optimum  performance. 

Progress  in  fulfilling  these  goals  is  summarized  briefly  in  the  following  paragraphs 
and  is  documented  in  detail  in  Appendix  II-A  (text  to  be  presented  at  the  AIAA  18th 
Aerospace  Sciences  Meeting,  January  1980). 

Tests  with  pure  base  burning  have  thoroughly  documented  the  encouraging 
base  pressure  rise  and  associated  performance  values  for  this  apparatus.  A  reliable 
firing  procedure  featuring  consumable  igniters  was  developed  and  repeatable 
measurements  were  obtained.  Wake  surveys  to  measure  pressure  and  temperature 
detail  were  completed  for  a  range  of  hydrogen  injection  rates.  Severe  deflection  of 
the  pitot-static  and  thermocouple  probes  during  these  tests  resulted  in  binding  and 
slippage  in  the  vertical  portion  of  the  traverse  program.  The  measurements, 
therefore,  are  questionable  for  determining  precise  spatial  pressure  or  temperature 
distributions.  The  temperature  measurements,  however,  are  valuable  for 
demonstration  of  the  lean-burning  nature  of  the  wake.  Temperatures  were  highest 
within  two  model  diameters  of  the  base  and  decreased  continuously  with  distance 
from  the  base.  In  all  cases,  wake  temperatures  remained  significantly  below  the 
flame  temperature  for  stoichiometric  burning  of  hydrogen.  Further  measurements 


with  the  present  probe  design  were  not  continued  so  that  "pure"  external  burning 
tests  could  be  pursued.  If  judged  valuable,  wake  detail  measurements  may  be 
resumed  in  the  coming  year  using  re-designed  probes. 

Attempts  to  achieve  pure  external  burning  have  not  yet  suceeded.  The 
original  radial  injection  configuration  using  six  orifices  recessed  in  an  annular 
channel  and  operating  subsonically  (M<0.1)  yielded  base  pressure  rises  and 
performance  values  similar  to  the  base  injection  tests.  Burning  was  obtained  in  the 
annular  channel  as  well  as  in  the  wake;  however,  no  combustion  in  the  free  stream 
was  evident.  Orifice  sizes  were  reduced  and  supersonic  injection  was  attempted  to 
increase  jet  penetration.  Ignition  with  jet  Mach  numbers  greater  than  about  2.1  was 
not  obtained  either  in  the  channel  or  in  the  wake.  With  injection  at  1.5  <M  <2.1, 
ignition  was  sustained  in  the  channel  but  not  in  the  wake  and  base  pressures  were 
not  significantly  elevated.  Further  tests  with  this  configuration  are  currently  in 
progress  to  explore  the  range  of  subsonic  and  supersonic  injection  conditions 
between  the  extremes  defined  by  the  original  tests. 

Investigation  of  combined  external/base  burning  was  effected  by 
superimposing  simulated  axisymmetric  and  asymmetric  external  burning  on  actual 
base  burning.  These  tests  revealed  that,  for  the  three  simulated  burning  fields 
applied,  the  resulting  base  pressure  is  very  nearly  the  simple  sum  of  the  individual 
base  and  external  contributions  when  acting  alone.  These  results  confirm  earlier 
tests  in  this  program  with  simulated  external  and  base  burning  showing  that  base 
pressures  can  be  elevated  substantially  above  free  stream  (i.e.,  base  thrust)  and  are 
consistent  with  those  recently  presented  by  Schadow  and  Chieze.^  Tests  exploring 
actual  combined  external/base  burning  are  being  conducted  concurrently  with  those 
attempting  to  achieve  "pure"  external  burning.  Radial  injection  at  decreasing  Mach 
numbers  is  combined  with  base  injection  through  sintered  plugs  of  varying  porosity 


to  maintain  an  approximately  equal  split.  The  success  of  these,  of  course,  will 
closely  parallel  the  success  of  tests  with  attempted  pure  external  burning. 

Preparation  for  the  coming  year's  proposed  base  burning  tests  of  diluted 
hydrogen  was  begun  in  the  current  year.  Design  calculations  were  completed  for 
purchase  of  flow  meters  to  obtain  the  desired  molecular  weight  ranges  with  the 
specified  diluent  gases.  Two  rotameters  each  with  two  alternate  floats  have  been 
purchased  and  delivered.  Gas  regulators  and  pressure  gages  have  been  purchased 
and  received.  All  other  plumbing  hardware  will  be  bought  at  the  time  of  assembly. 
Model  adjustmenets  for  these  tests  are  minor.  All  designs  are  complete  and 
fabrication  is  now  in  progress. 
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Abstract 

Investigations  continued  of  the  accumulation  processes  of  aluminum  on  the 
propellant  burning  surface  that  lead  to  formation  of  agglomerate  droplets,  the  size 
of  which  dominates  aluminum  combustion.  Using  dry-pressed  mixtures  of 
ammonium  perchlorate,  aluminum  and  carnauba  wax  powders,  the  effects  of 
particle  sizes,  mixture  ratios  and  pressure  were  studied  by  high  speed  photography. 
Results  were  correlated  with  a  scenario  describing  the  effects  of  propellant 
microstructure,  retention  processes  for  aiixninum  on  the  burning  surface,  response 
of  aluminum  powders  to  heating,  and  inf  lamation  process  for  sintered  aluminum. 

Quench  tests  at  atmospheric  pressure  on  the  combustion  plume  above  the 
burning  surface  showed  that  the  agglomerates  in  these  tests  burned  much  as 
aluminum  droplets  in  earlier  single  particle-controlled  atmosphere  tests.  An 
exception  resulted  in  tests  with  large  (95  ^m)  aluminum  particles  which  did  not 
ignite  completely.  An  apparatus  for  tests  at  elevated  pressure  was  designed  and 
partially  fabricated. 
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I.  Introduction 

(1.4) 

It  was  noted  in  earlier  progress  reports'  that  the  primary  unknowns  in 
combustion  behavior  of  aluminum  in  solid  rockets  arise  from  the  complicated 

4 

behavior  of  the  aluminum  on  the  propellant  burning  surface.  While  the  combustion 
occurs  primarily  after  the  aluminum  leaves  the  burning  surface,  the  nature  of  the 
burning  droplets  is  determined  by  surface  processes  involving  accumulation  and 
agglomeration  of  original  ingredient  particles.  The  present  studies  have  been 
aimed  at  determining  the  nature  of  these  processes,  the  nature  of  the  burning 
droplets,  and  the  effect  on  combustion  behavior. 

Initial  work  on  this  program  concentrated  on  the  behavior  of  the  propellant 
ingredients  during  heating,  including  formation  of  binder  melts  and  interaction  of 
particles  (particularly  aluminum).  These  studies  established  the  conditions 
necessary  for  concentration,  adhesion  and  sintering  of  particles  on  the  propellant 
surface.  Studies  were  then  made  of  the  evidence  of  these  processes  in  propellant^ 
like  systems  chosen  to  control  such  critical  conditions  as  the  presence  of  binder 
melts,  proximity  of  aluminum  particles  in  the  matrix  of  the  solid,  and  oxidizing 
character  of  the  surroundings  of  the  aluminum  particles.  These  studies  (on  samples 
made  from  AP/A1  mixtures)  suggested  that  molten  binder  was  not  essential  to 
retention  and  accumulation  of  aluminum  on  the  burning  surface  (although  it  was 
conducive  to  accumulation  under  some  conditions).  The  results  showed  that  the 
aluminum  was  resistant  to  ignition  in  the  AP  flame,  but  that  accumulating  particles 
would  sinter  together.  This  latter  process  was  favored  by  particle  size 
combinations  that  provided  proximity  of  aluminum  particles  to  each  other. 

Introduction  of  a  hydrocarbon  fuel  in  graduated  amounts  showed  that  ignition 
events  occurred  in  the  accumulating  aluminum,  with  surface-wise  spreading  of 


j 
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inflamation  sweeping  up  adjoining  aluminum  accumulations  to  give  burning 
agglomerates.  These  results,  and  tests  on  samples  with  thin  fuel  laminations  show 
that  diffusion  flamelets  associated  with  the  H.C.  fuel  vapors  act  to  pilot  the 
aluminum  ignition.  Tests  along  these  lines  have  been  continued  and  are  described 
below. 

At  the  beginning  of  this  report  period,  preliminary  tests  had  been  made  to 
further  evaluate  the  aluminum  behavior  after  leaving  the  burning  surface. 
Experiments  provided  samples  of  agglomerates  quenched  at  various  distances  from 
the  burning  surface.  The  limited  results  available  suggested  transformation  of 
accumulates  to  liquid  aluminum  droplets  near  the  surface,  followed  by  combustion 
to  burnout  without  any  singular  events. 

While  the  previously  reported  work  did  not  pretend  to  determine  details  of 
aluminum  behavior  at  the  molecular  or  subrmcron  level,  it  established  a  basis  for 
recognition  of  the  dominant  processes  governing  aluminum  behavior;  and  hence 
provided  some  rational  basis  for  modification  of  behavior.  More  recent  work  has 
sought  to  consolidate  the  basis  of  understanding  by  extending  previous  test  series 
and  by  development  of  a  mechanistic  model  embodying  the  understanding  to  date. 
In  addition,  observations  of  behavior  of  droplets  after  detachment  from  the  burning 


surface  have  been  continued. 


A.  Burning  Surface  Behavior:  AP/A1  Systems 

Various  studies  have  been  reported  on  combustion  of  solid  samples  made  by 
dry-pressing  mixtures  of  AP  and  A1  powders  There  is  a  tendency  for  such 

works  to  be  provincial  in  nature,  exploring  a  limited  range  of  experimental 
variables  and  making  a  broad  interpretation  based  on  the  rather  limited  range  of 
observational  data  (e.g.,  argument  about  aluminum  behavior  based  on 
determinations  of  sample  burning  rate).  In  the  present  studies,  more  direct 
observation  of  the  samples  was  made  using  high  speed  cinemicrophotography.  Test 
variables  were  chosen  to  reflect  rocket  motor-like  conditions  (within  the  limits  of 
the  window-bomb  type  of  experiment).  Observations  were  also  made  of  samples 
quenched  by  various  methods.  State  and  behavior  of  the  aluminum  itself  was 
addressed  in  these  observations.  Work  during  this  reporting  period  on  the  AP/A1 
system  included  extension  of  the  previous  test  series  and  more  thorough 
examination  of  results.  Table  III-l  summarizes  the  complete  series  and  test 
conditions.  In  many  cases  the  tests  are  duplicates,  run  to  obtain  more  decisive 
information  or  better  motion  pictures. 

The  most  productive  tests  were  the  motion  pictures.  These  were  viewed 
independently  by  several  observers,  who  evaluated  the  behavior  relative  to  the 
different  aspects  of  behavior  listed  in  Table  III-2.  After  recording  these 
evaluations,  the  observers  compared  results  and  agreed  on  the  following  observa¬ 
tions. 

1.  Aluminum  accumulates,  concentrates  and  sinters  on  the  burning  surface 
before  detaching. 

2.  Samples  containing  AP  and  A1  of  comparable  particle  size  showed  only 
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Table  II1-2 

Topics  to  Evaluate  From  Movies 
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1.  Degree  of  surface  accumulation  (s.a.) 

2.  Connectedness  of  s.a. 

3.  Looseness  of  attachment  of  connected  s.a.  (c.s.a.)  (i.e.,  apparent  ri¬ 
gidity  vs.  freedom  of  movement  while  attached) 

crust?  filigree? 

4.  To  what  extent  do  aluminum  particles  or  accumulates  ignite  in  the 
field  of  view? 

5.  Is  the  ignition  primarily  before,  during  or  after  detachment 

6.  Are  most  of  the  particles  leaving  the  surface 

original  A1  particles 

accumulates  of  Al  particles  (typically  how  many?) 
ignited  agglomerates 

7.  Describe  size  and  size  distribution  of  burning  droplets. 

8.  How  vigorous  is  the  Al  combustion?  "Vigorous"  is 

many  small  bright  droplets 
overall  plume  bright 

smoke  density  in  plume  builds  up  close  to  surface  (don't  count 
recirculating  smoke,  etc.) 

9.  Is  flux  of  Al  from  the  burning  surface  macroscopically  continuous; 
spasmodic  (in  time  or  location);  periodic;  surface  propagative  (e.g., 
in  accumulate  layer)? 


limited  accumulation  and  sintering.  Coarse  AP  with  fine  A1  showed  massive 
accumulation  and  sintering. 

3.  The  environment  of  the  AP  flame  is  not  effective  in  igniting  the  aluminum 
in  the  condition  present  on  the  burning  surface.  In  the  case  of  coarse  A1  particles, 
most  left  the  surface  without  igniting  (i.e.,  without  visible  flame  envelopes).  The 
sintered  accumulate  obtained  with  fine  particles  did  not  ignite  readily,  and  when  it 
did  ignite,  the  conditions  suggested  that  break-up  of  the  accumulate  played  a 
causative  role  (i.e.,  onset  of  mobility  and  detachment  seemed  to  precede  ignition). 

4.  When  ignition  occurred  in  the  sintered  accumulate,  it  tended  to  spread 
through  the  accumulate.  In  the  process,  the  accumulate  layer  coalesced  into  a 
number  of  large  burning  droplets  (referred  to  here  as  "agglomerates"),  in  a 
progressive  process  involving  also  progressive  detachment  from  the  surface. 

5.  Aluminum  that  experienced  detachment  without  ignition  usually  did  not 
ignite  later.  In  the  case  of  sintered  fine  AJ,  sintered  structures  broke  up  some 
during  detachment,  and  those  fragments  that  were  already  ignited  somewhere  in 
the  fragment  became  fully  inflamed. 

6.  The  behavior  described  above  was  affected  only  moderately  by  pressure 
over  the  pressure  range  tested  (4.1  to  6.9  MPa,  600-1000  psi). 

7.  Fine  ingredient  particle  size  (10/3)  led  to  sporadic  inflamation  and 
agglomeration  at  the  surface;  inflamation  and  combustion  of  the  A1  was  then 
extremely  vigorous. 

In  terms  of  mechanistic  arguments  relative  to  propellant  combustion,  the 
observations  1-7  (combined  with  results  reported  earlier  on  behavior  of  A1  particles 
during  heating)  imply  the  following! 

*  The  numbered  mechanistic  arguments  pertain  to  the  correspondingly 

numbered  observations  of  the  preceding  list. 


1.  Whatever  the  state  of  the  aluminum  particles  and  the  AP  surface  when  the 
AP  particles  emerge  on  the  burning  surface,  there  is  some  adhesive  process  that 
prevents  immediate  detachment  of  the  particles.  Under  suitable  conditions  (e.g., 
AP  size  large  compared  to  A1  size),  retention  is  long  enough  to  result  in  emergence 
of  more  A1  particles,  concentration  of  particles  and  sintering.  Retention  may 
initially  involve  the  liquid  state  of  the  AP  surface,  and  conditions  at  the  surface 
suggest  that  sintering  is  similar  to  that  observed  slightly  above  the  aluminum 
melting  point  in  heating  tests  on  A1  powders  in  oxidizing  atmospheres. 

2.  Formation  of  aggregates  of  metal  particles  depends  on  both  the 
mechanism  and  the  opportunity.  Work  to  date  has  focused  on  mechanisms  by  which 
aluminum  particles  stick  to  each  other  (sintering,  coalescence),  and  the  necessity 
and  means  for  adhesion  to  the  burning  surface  to  permit  particle  concentration 
(opportunity  for  particle  interaction).  However  the  opportunity  is  governed  also  by 
how  long  a  particle  must  "wait"  on  the  surface  before  it  is  joined  by  other  particles, 
and  by  the  nature  of  its  environment  during  the  wait.  Thus  opportunity  is 
dependent  on  microstructure  of  the  sample,  including  particularly  the  grouping  of 
aluminum  in  pockets  and  connective  arrays  between  pockets.  When  the  number  of 
oxidizer  and  aluminum  particles  are  equal,  the  opportunity  for  aluminum  particle 
interaction  is  low,  the  wait  for  interaction  long,  and  the  retention  environment 
poor,  so  that  sintering  is  minimal.  When  the  aluminum  particle  size  is  small,  the 
particles  group  into  packing  voids  among  AP  particles,  and  may  even  form 
connective  arrays  among  the  voids.  This  provides  the  opportunity  for  accumulation 
and  sintering  of  groups  of  particles  and  interconnecting  bridges  between  groups, 
giving  the  sintered  layers  observed  in  tests. 

3.  Mechanistically,  the  reluctance  of  aluminum  to  ignite  in  most  of  these 
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tests  indicates  that  the  oxide  coating  on  the  particles  is  resistive  to  breakdown  in 
the  chemical  and  thermal  environment  of  the  AP  flame.  This  is  analogous  to 
results  of  hot  plate  experiments  where  agglomeration  of  aluminum  powders 
occurred  at  or  near  to  the  aluminum  melting  point  in  an  Argon  atmosphere,  but  not 
in  an  oxygen  enriched  one.  It  appears  that  the  oxide  "skin"  is  reinforced  before 
ignition,  even  in  the  AP  flame,  thus  delaying  ignition.  This  exceptional  behavior  of 
aluminum  (and  some  other  metals)  is  unique  among  fuels:  oxidation  temporarily 
inhibits  combustion.  In  any  case,  the  protective  property  of  the  oxide  skin  in 
conventional  propellants  contributes  to  the  opportunity  for  sintering  and 
agglomeration,  by  delaying  ignition,  even  in  the  AP  flame. 

4.  From  observations  3  and  4  above,  it  is  concluded  that  sintered 
accumulates  consist  of  interconnected  oxide  shells  containing  molten  aluminum 
that  is  exposed  when  the  sintered  structure  is  broken.  This  can  lead  to  local 
inflamation  of  aluminum,  with  corresponding  local  temperature  rise.  Given  the 
relatively  low  thermal  conductivity  of  a  sintered  accumulate  (compared  to 
aluminum),  the  local  temperature  rise  is  sufficient  to  melt  the  surrounding  oxide, 
permitting  propagation  of  the  ignition  front  through  the  accumulate.  The 
aluminum  is  drawn  by  surface  tension  into  one  or  more  burning  droplets. 

5.  The  interpretation  in  4  is  further  supported  by  the  observation  in  5  that 
ignition  rarely  occurs  in  single  or  sintered  particles  if  they  do  not  ignite  during  the 
mechanical  stresses  of  detachment  from  the  burning  surface.  However,  this 
support  must  be  tempered  by  recognition  that  the  tests  were  made  with  small 
samples,  that  slow  reaction  probably  continues  without  flame  type  combustion,  and 
that  other  investigators  have  obtained  intense  combustion  using  very  fine  AP  and 
A1  particle  size. 


/() 
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6.  The  insensitivity  of  behavior  to  pressure  may  be  in  part  related  to  the 
limited  pressure  range  tested  (4.1  to  6.9  MPa).  Lower  pressures  were  not  tested 
because  the  samples  would  not  support  combustion.  This  is  a  significant  point 
because  other  investigators  have  gone  to  much  lower  pressures,  but  used  very  fine 
AP  and  Al.  Their  observations  regarding  the  role  of  condensed  phase  reactions  in 
determining  burning  rate  are  apparently  not  relevant  to  samples  made  from 
propellant-size  ingredients.  Such  reactions  are  involved  in  sintering  and  in  delay  of 
aluminum  ignition,  and  are  not  conspicuously  pressure-dependent. 

7.  Behavior  with  very  fine  ingredients  has  been  studied  by  others  as  well,  and 
collected  results  suggest  that  substantial  surface  reaction  of  Al  particles  takes 
place  at  the  propellant  burning  surface  due  to  the  high  specific  surface  of  the 
reactants.  This  may  involve  burning  off  small  areas  of  the  heated  AP/A1  surface 
layer  rather  than  outright  aluminum  accumulation.  In  that  case  the  arguments  for 
non-accumulation  due  to  Al  particle  isolation  in  the  microstructure  may  be 
irrelevant  to  the  tests  with  10/5  mixtures. 

B.  Burning  Surface  Behavior  in  AP/A1  Systems  With  Modified  Aluminum 

In  view  of  the  obvious  important  role  of  the  oxide  skin  in  controlling  the  onset 

of  sintering,  agglomeration  and  ignition  of  aluminum  particles,  it  is  reasonable  to 

(9) 

seek  beneficial  modification  of  the  oxide.  A  method  explored  by  Krauetle  was  to 
enhance  the  oxide  by  further  oxidation  by  holding  powders  at  elevated  temperature 
in  oxidizing  atmospheres.  This  method  was  called  "pre-oxidation",  and  was 
conducted  at  temperatures  below  the  aluminum  melting  point. 

A  method  explored  earlier  in  the  present  project  was  called  "pre-stretching" 
the  oxide  by  heating  particles  through  the  aluminum  melting  point.  The  oxide  skin 
deforms  to  accomodate  the  increased  volume  of  the  aluminum,  probably  by  both 
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inelastic  stretching  and  cracking.  In  the  presence  of  a  low  concentration  of 
oxygen,  the  cracked  areas  will  close  rapidly  by  further  oxidation.  Upon  cooling,  the 
particles  will  shrink,  the  oxide  skin  will  wrinkle  or  exhibit  depressions,  but  the 
oxide  surface  area  will  remain  sufficient  to  enclose  the  aluminum  when  the  particle 
melts  in  the  combustion  zone.  This  argument  was  developed  from  growing 
understanding  of  preignition  behavior  of  aluminum  powder,  and  was  evaluated 
earlier  in  the  project  using  the  hot  stage  microscope  to  produce  and  test  the  pre¬ 
stretched  oxide  particles.  In  those  tests  the  tendency  of  aluminum  powders  to 
sinter  and  agglomerate  when  heated  was  sharply  reduced  by  pre-stretching  the 
oxide. 

In  the  present  report  period,  aluminum  with  pre-stretched  oxide  was  produced 
in  greater  quantity  by  heating  the  powder  in  an  open  quartz  tube  to  700°C,  using  a 
tube  furnace  flushed  with  a  nitrogen  flow  (with  some  entrained  air).  Evaluation 
was  then  made  oy  preparation  and  combustion  of  dry-pressed  samples  of  85%,  100 
um  AP  and  15%  Alcoa  123  Al.  Similar  tests  were  also  run  on  samples  with  untreated 
Alcoa  123  Al,  and  samples  with  pre-oxidized  Alcoa  123A1  (provided  by  Karl  Kraeutl** 
of  Naval  Weapoas  Center).  Tests  were  run  at  6.9  MPa  (1000  psi),  and  observations 
were  made  by  combustion  photography  (Table  III-3). 

Tests  on  the  samples  with  untreated  aluminum  were  similar  to  earlier  ones 
(Table  III-l)  with  100  um  AP  and  H-30  aluminum,  exhibiting  massive  accumulation 
and  sintering  of  aluminum  on  the  burning  surface,  with  ignition  occurring  only 
during  break-up  of  detaching  accumulate  layers.  Very  large  agglomerates  formed. 
Results  with  the  pre-oxidized  and  with  the  pre-stretched  aluminum  were  alike. 
Only  small  accumulates  were  evident,  with  more  or  less  continual  detachment  of 
fragments.  Aluminum  ignition  was  only  occasional.  This  result  supports  the 
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Evaluation  of  Aluminum  Pretreatment  Table  Identifies  Combustion  Photography 
Tests.  All  samples  from  dry-pressed  mixtures  of  85%  100  pm  AP  and  15% 
Alcoa  123  aluminum  tested  at  6.9  MPa. 


Type  of 

Aluminum 

Treatment 

None 

Pre-oxidized 

Pre-stretched 

63 

66, 

65,  67 

Film  # 

64,  70 

69,  72 

71 

Test 

Fair  Quality 

Fair  Quality 

Fair  Quality 

Outcome 

(smoke) 

(smoke) 

(smoke) 

mechanistic  argument  that  led  to  ''pre-stretching*'  experiments,  and  suggests  a 
means  of  controlling  accumulate  size,  using  a  modification  of  aluminum  powder 
that  is  economically  viable  in  production,  possibly  by  simply  changing  process 
control  variables  in  the  original  powder  manufacture. 

C.  Burning  Surface  Behavior:  AP/Al/Wax  Systems 

In  view  of  the  reluctance  of  aluminum  to  ignite  in  the  AP  flame,  it  was 
concluded  that  the  diffusion  flamelets  from  the  oxidizer-binder  vapors  played  a 
role  in  aluminum  ignition  in  propellant  flames  (a  postulate  advanced  earlier  in 
connection  with  tests  on  "sandwiches"  with  aluminized  binders®  and  other  tests 
with  aluminized  AP It  has  often  been  suggested  that  aluminum  ignition  would 
not  be  effectively  accomplished  unless  the  hydrocarbon  fuel-AP  (HC-AP)  mixture 


ratio  was  suitable  for  an  equilibrium  combustion  temperature  above  the  melting 
point  of  aluminum  oxide,  thereby  assuring  breakdown  of  ail  protective  oxide. 
However,  our  observations  of  the  local,  propagative  nature  of  the  ignition  of 
accumulate  layers  suggest  that  satisfactory  ignition  would  occur  if  a  reasonable 
distribution  of  local  HC-AP  diffusion  flames  occurred,  even  though  that  might 
correspond  to  a  relatively  small  mass  fraction  of  HC  fuel. 

Since  the  accumulation  process  (and  hence  agglomerate  size)  tends  to  be 
limited  by  local  ignition  events,  it  seems  important  to  understand  what  local 
diffusion  flamelet  situations  are  sufficient  for  ignition.  This  was  explored  by 
adding  a  hydrocarbon  (carnauba)  wax  to  the  AP/Ai  mixture  for  dry-pressed  samples. 
Combustion  photography  tests  were  run  that  were  similar  to  those  for  AP/AI 
samples.  The  complete  series  of  tests  is  identified  in  Table  III-4,  and  the  test 
results  were  reviewed  in  the  same  manner  as  the  AP/AI  tests.  Observations  were 
as  follows: 

1.  Addition  of  2%  HC  wax  reduced  surface  accumulation  to  a  mininum,  giving 
a  steady  flow  of  mostly  nonburning  A1  from  the  surface. 

2.  4%  wax  resulted  in  increased  surface  accumulation  of  A1  (much  less  than 
for  0%  wax,  comparable  to  propellants),  with  sporadic  detachment  and  considerable 
inflamation. 

3.  8%  wax  resulted  in  formation  of  burning  agglomerates  on  the  surface, 
which  immediately  detached  and  burned  vigorously  in  the  outward  flow. 

4.  In  all  cases,  increased  pressure  led  to  smaller  accumulates  and  more 
extensive  ignition  of  aluminum  on  or  near  the  surface. 


Summary  of  Tests  on  Dry -Pressed  Ap  /  A1  /  Wax  Mixtures* 
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D.  Theory  o f  Accumulation-Agglomeration 

The  processes  involved  in  the  behavior  of  aluminum  on  the  burning  surface 
are  clearly  rather  intricate,  and  may  never  be  well  described  by  a  quantitative 
theory.  However  it  is  important  to  develop  an  approximate  theory,  in  order  to 
pursue  intelligent  design  of  experiments  and  to  fully  exploit  results  of  past 
experiments  in  both  future  experiments  and  practical  applications.  Certain  steps  in 
the  aluminum  behavior  have  been  defined  in  order  to  assure  that  relevant  processes 
are  considered,  and  the  relation  of  these  steps  are  identified  in  Figure  1II-1  and 
discussed  below.  Items  1  and  2  are  initial  conditions,  while  the  others  are 
processes. 

1.  First  are  the  many  questions  about  the  nature  of  the  individual  propellant 
ingredients,  such  as  the  condition  of  the  oxide  coating  on  the  aluminum  particles. 

2.  Second  is  the  microstructure  of  the  propellant,  which  reflects  such 
variables  as  particle  size  or  specific  surface,  proximity  of  A1  particles  to  each 
other,  and  nature  of  the  material  near  the  aluminum. 

3.  Third,  as  the  burning  surface  approaches  one  or  a  group  of  aluminum 
particles,  the  temperature  rises  until  the  surrounding  binder  is  softened  and 
gasifying.  In  effect,  the  aluminum  particle  is  wetted  by  melting  binder,  and 
scattered  particles  may  be  concentrated  in  the  •'melt"  as  binder  is  removed  by 
vaporization.  In  a  heterogeneous  propellant,  this  behavior  is  statistically 
distributed. 

4.  Aluminum  particles  must  each  be  reached  by  the  regressing  surface 
between  condensed  and  gas  phase.  At  that  point  they  evidently  are  retained  on  the 
surface,  by  adhesion  to  the  melt  and/or  underlying  Ai  particles,  leading  to 
continued  accumulation.  Ignition  may  be  delayed  by  a  continuing  fuel-rich 
environment,  and  by  protectiveness  of  the  oxide  on  the  particles  at  temperatures 
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present  at  the  burning  surface. 

5.  Increasing  interaction  among  the  adjoining  particles  leads  to  sintering  by 
either  carbonaceous  material  or  direct  particle  interaction  (oxidation),  or  both  in 
that  order.  Temperatures  at  this  point  may  exceed  the  aluminum  melting  point. 
The  microstructure  of  the  propellant  (especially  spatial  distribution  of  aluminum) 
plays  a  key  role  in  determining  the  opportunity  for  particles  to  encounter  each 
other  in  or  on  the  burning  surface,  and  in  determining  the  extent  and  nature  of 
sintering.  Exposure  to  an  increasingly  oxidizing  atmosphere  can  cause  sintering 
without  causing  ignition. 

6.  It  is  helpful  to  consider  the  formation  of  sintered  accumulates  in  two  or 
three  ways.  In  its  simplest  form,  very  small  numbers  of  particles  may  contact  each 
other  before  detachment  from  the  surface.  A  more  deterministic  process  is  the 
accumulation,  sintering  and  detachment  of  all  particles  in  each  packing  void  in  the 
AP  array  (assemblages  of  small  aluminum  particles,  created  during  preparation  of 
the  propellant).  This  description  of  accumulate  size  distribution  is  called  the 
pocket  model.  A  more  complete  model  would  have  to  consider  formation  of  bridges 
between  pocket  accumulates,  leading  to  filigrees,  flakes  or  layers.  Bridging  is 
related  to  the  nature  and  aluminum  content  of  the  passages  between  voids  in  the 
AP  packing  array. 

7.  Agglomeration  and  ignition  are  usually  closely  coupled  events,  both  being 
governed  by  breakdown  of  oxide  and  expo  sire  of  aluminum.  To  describe  the 
occurrence  on  a  propellant  surface,  one  must  first  describe  the  development  of 
accumulates  and  exposure  to  diffusion  flamelets,  and  then  consider  how  an  ignition 
site  spreads  through  the  accumulate  to  form  one  or  more  agglomerates,  and 
consider  what  limits  the  extent  of  this  surface-wise  inflamatipn.  In  an  ordinary 
propellant,  a  proper  statement  of  this  ignition-inf  lamation  problem  would  normally 
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have  to  be  made  in  statistical  terms  related  to  propellant  microstructure.  To  get 
started,  however,  one  considers  the  inflamation  process  through  a  uniform  sintered 
structure,  and  then  applies  the  results  to  a  non-uniform  one  having  the 
characteristics  of  a  propellant  surface. 

8.  The  extent  of  the  accumulation-sintering-agglomeration  process  is 
dependent  always  on  opportunity,  and  is  statistically  variable  in  more  ways  than  the 
surface-wise  inflamation.  Aluminum  detaches  from  the  burning  surface,  at  which 
point  the  opportunity  for  further  sintering  or  agglomeration  is  drastically  reduced. 
Detachment  results  from  the  interplay  of  restraining  and  removing  forces  at  the 
surface,  which  change  progressively  during  accumulation  and  differ  greatly  from 
point  to  point  on  the  surface,  from  one  propellant  to  the  next,  and  from  one  flow 
environment  to  the  next.  Description  of  these  forces  has  been  attempted  for  only 
the  most  simplistic  situations  (e.g.,  spheres  in  an  evaporating  viscous  layer).  At 
present  it  seems  more  important  to  identify  the  nature  of  an  accumulate,  the 
nature  of  its  attachments  to  a  real  propellant  surface,  the  way  in  which  the 
ignition-agglomeration  process  changes  the  interplay  of  forces,  and  what  conditions 
precipitate  ignition  (in  most  AP-HC  binder  propellants,  detachment  is  precipitated 
by  ignition,  and  growth  of  the  agglomerate  is  terminated  by  completion  of 
detachment). 

A  comprehensive  theory  of  accumulation-agglomeration  would  have  to 
encompass  the  concepts  1-8.  The  experimental  work  has  been  aimed  at  identifying 
the  most  important  steps  and  clarifying  them.  Analytical-computational  studies  of 
flame  spreading  from  ignition  sites  have  also  been  made  in  the  present  program. 
These  studies  were  aimed  at  describing  periodic  accumulation-ignition-inflamation 
observed  in  some  AP/A1  combustion  tests.  A  second  study  has  examined  the  void 
structure  resulting  from  packing  particles  (i.e.,  the  microstructure  of  the  volume 


filled  by  binder  and  aluminum  in  an  AP  propellant).  The  goal  of  this  latter  study  is 


to  determine  appropriate  statistical  variables  to  characterize  the  distribution  of 


fine  aluminum  particles  in  the  propellant  microstructure.  These  two  studies  are  in 


early  phases. 


E.  Quench  Studies  of  the  Combustion  Plume 


In  a  rocket  motor,  the  majority  of  the  aluminum  in  the  propellant  burns  as 


droplets  while  moving  out  with  the  flow  from  the  propellant  surface.  Because  of 


the  agglomeration  process,  these  droplets  are  large  enough  to  require  appreciable 


time  to  burn.  In  extreme  cases,  droplets  may  not  be  completely  burned  at  the 


nozzle  entrance,  in  which  case  loss  of  efficiency  occurs.  In  all  cases,  there  exists  a 


significant  volume  containing  burning  droplets,  a  volume  that  also  participates  in 


erosive  burning  processes  and  in  oscillatory  combustion-damping  phenomena.  The 


effects  of  this  participation  have  not  been  addressed  for  a  variety  of  reasons,  the 


most  basic  being  ignorance  of  the  characteristics  of  the  droplet  cloud.  This  subject 


was  addressed  in  another  project^  in  which  some  success  was  achieved  in  quench- 


collection  of  aluminum  droplets  at  various  distances  from  the  burning  surface 


(atmospheric  pressure). 


The  scenario  for  agglomerate  combustion  has  been  described  before,  ’  and 


(10,11) 


the  object  of  the  tests  was  to  obtain  more  decisive  evidence,  and  explore  the 


limitations  of  the  scenario.  It  is  widely  (but  not  uniformly)  believed  that  aluminum 


agglomerates  burn  like  aluminum  droplets,  with  oxidation  in  a  detached  flame  and 


on  the  droplet  surface.  Products  are  predominantly  AlpOj  smoke  from  the 
detached  flame,  with  perhaps  20%  of  the  AljO^  forming  on  the  droplet  surface  and 
yielding  one  or  more  relatively  large  "residual"  oxide  droplets  after  aluminum 


burnout.  In  this  scenario,  a  sample  quenched  from  the  plume  should  contain  smoke 


oxide  to  pm  diameter),  residual  oxide  (typically  10-60  pm  diameter),  and  aluminum 
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particles  that  are  partly  burned  and  have  oxide  lobes  that  are  the  source  of 
"residual"  oxide  particles  upon  droplet  burnout.  The  size  of  the  partially  burned 
aluminum  droplets  depends  on  the  degree  of  agglomeration  on  the  burning  surface 
(i.e.,  original  agglomerate  size)  and  the  distance  out  to  the  quench  site. 

The  above  scenario  for  agglomerate  combustion  is  based  on  contributions 
from  many  sources,  very  few  of  which  were  known  to  be  directly  applicable  to 
propellant  combustion  (which  is  very  difficult  to  observe  unambiguously). 
Outstanding  practical  questions  to  be  addressed  in  the  present  and  future  work  are: 

1.  Nature  of  the  original  agglomerate,  and  relevance  of  single  particle 
combustion  studies  on  pure  aluminum  (this  includes  the  question  of  how  much  oxide 
is  already  present  on  the  agglomerate  by  the  time  it  is  fully  inflamed  ("pre-ignition 
oxide")). 

2.  Amount  of  additional  oxide  accumulation  on  the  agglomerate  during 
burning. 

3.  Effect  of  surface  oxide  on  the  burning  history  and  burning  time. 

4.  Extent  of  subdivision  of  surface  oxide,  either  by  ejection  during  burning  or 
division  during  droplet  burnout. 

5.  Effect  of  a  convective  flow  around  the  burning  agglomerate. 

6.  Responsiveness  of  droplet  burning  to  oscillatory  flow. 

7.  Dependence  of  1-6  on  pressure. 

8.  Dependence  of  1-7  on  chemistry  and  temperature  of  the  atmosphere,  and 
radiation  field  in  the  rocket  motor  environment. 

In  the  present  studies,  the  quench-collection  procedure  in  Ref.  10  was  used  to 
sample  the  combustion  plume  of  a  variety  of  model  propellant  formulations  (dry- 
pressed),  and  the  collected  samples  were  studied  to  clarify  the  details  of  the 
combustion  plume.  Tests  were  run  at  atmospheric  pressure,  and  construction  was 
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started  on  a  test  chamber  for  testing  at  elevated  pressure.  Table  III- 5  shows  a 
summary  of  the  tests  conducted.  Variables  were: 

Aluminum  Particle  Size 

Alcoa  123  ~  25  Mm 

H-95  ~  95  Mm 

Ammonium  Perchlorate  Particle  Size 
Fine  ~  5M>m 

Coarse  ~  100  Mm 

%  A1 

2,  4,  8, 12, 15  f 

AP/Wax  ratio  (HC  fuel  content) 

88/10,  88/3 

Quenching  Distance  (cm) 

5.3, 10.6, 15.9,  31.8 

"Q"  entries  in  the  table  indicate  conditions  for  which  quench  tests  were  run,  "M" 
designates  motion  pictures. 

The  results  of  the  combustion  experiments  can  be  summarized  as  follows: 

1.  Quench  samples  showed  the  smoke  oxide,  residual  oxide  and  agglomerates 
described  in  the  classical  scenarios  above,  in  all  tests  where  agglomeration  was 
possible.  Coarse  aluminum  yielded  also  "large”  particles  of  irregular  shape  that 
were  judged  to  be  original  particles  that  did  not  spheroidize  but  burned  slowly  to 
yield  substantial  surface  oxide.  In  tests  showing  only  these  irregular  aluminum 
particles,  no  other  oxide  was  evident. 

2.  Agglomerate  particles  obtained  in  this  experiment  appeared  black  with 
dull,  slightly  grainy  surfaces.  They  decreased  in  size  with  quenching  distance,  and 


Table  HI-5d 

Effect  of  Quenching  Distances 
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usually  develop  oxide  lobes  at  increasing  quench  distance.  The  agglomerates  were 
spherical  except  lor  the  lobes  and  some  other  anomalies  attributed  to  the  quench 
process  and  subsequent  HC1  acid  attack.  in  samples  with  low  aluminum 
concentration  (2  and  4%)  agglomerates  are  very  small.  In  all  tests,  agglomerates 
were  largely  consumed  at  7.6  cm.  Similar  spherical  particles  occurred  in  some 
tests  involving  coarse  aluminum,  with  similar  evolution  of  form  with  distance. 
These  "agglomerates"  are  evidently  the  classical  burning  droplets  with  oxide  lobes. 
Their  detailed  structure  was  not  determined,  but  they  dissolved  in  10%  HC1  in 
water  leaving  the  oxide  lobes. 

3.  Residual  oxide  spheres  (translucent,  10-60  urn  diameter)  were  produced 
only  in  combustion  plumes  having  the  type  of  spherical  particles  described  in  2 
above,  presumed  to  be  fully  burning  droplets.  About  one  in  five  residual  oxide 
spheres  is  shiny  black  in  character;  these  are  believed  to  have  inclusions,  possibly 
of  aluminum.  The  oxide  spheres  did  not  dissolve  in  acid.  The  number  and  size  of 
the  spheres  increases  with  distance  from  the  burning  surface,  and  follows  the  same 
trend  as  "agglomerate"  size  with  increase  in  %  A1  in  the  sample. 

4.  Oxide  smoke  was  conspicuously  present  in  those  tests  that  had 
"agglomerate"  spheres,  was  minimal  or  absent  in  tests  yielding  exclusively  non- 
spheroidized  aluminum.  No  attempt  was  made  to  study  the  smoke  oxide  because  it 
was  only  partially  recovered  in  the  tests,  and  tended  to  obscure  the  rest  of  the 
particles.  The  smoke  was  mostly  removed  with  the  HC1  in  the  washing  operation. 
The  presence  of  smoke  exclusively  with  tests  having  "agglomerate"  spheres  is 
evidence  that  the  non-spherical  aluminum  particles  have  burned  without  detached 
flame  envelope. 

5.  In  general,  the  observations  conform  with  the  results  of  the  high  pressure 
combustion  photography  tests  relative  to  agglomeration  of  aluminum  on  the  burning 


surface.  Coarse  aluminum  particles  leave  the  burning  surface  individually,  fine 
particles  agglomerate  if  present  in  higher  concentrations.  Fine  AP,  with  fine  A1  in 
lower  concentrations  gives  very  vigorous  A1  combustion. 

6.  Photographic  observations  revealed  no  fragmentation  events,  but  technical 
difficulties  have  continued  to  prevent  observation  of  the  burnout  phase  of 
agglomerates  unambiguously. 

In  general,  the  quench  tests  supported  the  scenario  for  agglomerate 
combustion,  except  in  the  case  of  coarse  A1  particles,  which  usually  did  not  get  hot 
enough  to  melt  the  oxide  and  give  vigorous  burning  at  this  pressure. 
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Figure  III-l.  Progression  of  aluminum  behavior  on  the  propellant  surface. 
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Abstract 

Facility  modification  was  completed  and  extensive  testing  was 
performed  on  an  apparatus  to  investigate  turbulence  induced  pressure 
fluctuations  in  a  rocket-like  cavity.  To  overcome  prior  deficiencies  in 
theory,  a  new  theory  was  constructed  for  the  generation  of  pressure 
disturbances  by  the  turbulence.  Undeniable  separation  has  been  achieved 
between  propagationai  and  local  pressure  fluctuations,  and  the  agreement 
between  theory  and  experiment  are  adequate.  Measurements  include 
microphone  and  hot  film  anemometry,  and  time  series  analysis  is  applied  to 
the  results. 

I.  Introduction 

The  background  pressure  fluctuation  level  in  a  rocket  motor  (liquid  or 
solid)  is  generally  of  the  order  of  1-2%,  even  under  nominally  stable 
operating  conditions.  These  fluctuations  are  random  noise,  modified  by  the 
chamber  acoustics  to  favor  frequencies  near  acoustic  resonant  mode 
frequencies.  This  background  noise  level  is  the  source  of  several  real  or 
potential  problem  areas.  It  is  a)  a  source  of  motor  vibration,  b)  a  possible 
disturbance  for  determining  the  location  of  transition  from  laminar  to 
turbulent  flow  in  the  chamber,  and  c)  a  possible  disturbance  to  trigger 
nonlinear  instability. 

Motor  vibration  interacts  in  a  detrimental  way  with  guidance 

accuracy  and  missile  weight.  Somewhat  recent  problems  with  Poseidon  Stage 

II  attest  to  this  fact/*)  More  recently,  troubles  have  occurred  with  the 
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space  shuttle  solid  rocket  booster,  However,  a  cause  of  the  motor 


vibration,  internal-pressure  fluctuations,  has  not  been  investigated,  nor  is 


the  background  fluctuation  level  usually  measured.  The  reason  for  a  lack  of 
information  is  that  pressure  instrumentation  is  usually  set  up  with  a 
sufficiently  low  gain  to  capture  large  amplitude  instabilities  -  not  noise.  The 
relatively  intense  acoustic  background  noise  is  the  source  of  the 
disturbances  which  may  be  amplified  if  the  motor  is  unstable.  Moreover,  it  is 
possible  that  even  if  the  motor  is  not  unstable  in  a  linear  sense  it  may  be 
driven  unstable  by  the  background  noise.  For  these  reasons  it  is  desirable  to 
characterize  this  noise  and  its  cause. 

Prior  work  on  this  and  another  contract  have  pointed  to  the  natural 
turbulence  in  the  rocket  cavity  gases  as  the  culprit  in  generation  of 
unwanted  pressure  fluctuations.  The  purpose  of  this  program  is  to  quantify 
the  relation  between  the  turbulence  and  the  generated  pressure  field. 

11.  Program  Progress 

The  details  of  the  past  year's  progress  are  located  in  Ref.  (3), 
attached  as  App.  IV- A  of  this  report.  The  tasks  completed  during  the  past 
year  are  as  follows: 

a)  Modification  of  the  facility  to  improve  signal  to  noise  ratio 
and  characterization  of  the  signal  into  propagational  sound  and 
local  noise 

b)  Extensive  testing  using  microphone  and  hot  film  anemometry 
measurements.  Time  series  analysis  has  been  applied  to  the 
results 


c)  Generation  of  a  new  theory  of  the  aeroacoustics  of  the  device 
which  explains  the  experimental  results. 
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At  the  close  of  the  contract  year  modifications  were  being  made  to 
the  apparatus  to  allow  side  wall  blowing,  to  simulate  combustion  gas  inflow 
from  the  sides. 
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Thin  paper  dcsctihcs  the  development  and  application  of  the  impedance  tube  technique 
in  the  measurement  of  localized  combustion  process  admittances.  The  toller  «ie  nested  for 
combustor  suability  analyses  The  cxpciimcntal  setup  consists  of  a  lube  with  the  com¬ 
bustion  process  to  be  investigated  at  one  end  and  an  acoustic  driver  at  the  other  end.  The 
latter  it  used  to  excite  a  standing  wave  in  the  tube  whose  struetute  depends  upon  its 
interaction  «‘<h  the  combustion  jcrocctt.  The  needed  combustion  process  admittance  it 
determined  ftom  acoustic  pressure  measutementt  taken  at  several  l«>c» lions  along  the  tube. 
The  paper  discusses  the  theoretical  foundation  of  the  ct  pci  intent  and  the  developed  data 
reduction  procedure  which  accounts  for  the  presence  of  steady  slate  temperature  gradient, 
mean  now  gradient  and  gas-phase  losses  in  the  tube.  Results  showing  the  applicability  of 
the  developed  experimental  technique  are  provided. 


».  INTRODUCTION 

This  paper  describes  the  adaptation  of  the  impedance  lube  technique  in  the  measurement 
of  the  admittances  of  localized  combustion  processes.  The  discussed  measurement  tech¬ 
nique  has  been  developed  as  part  of  an  ATOSR  teseutch  grant  that  was  concerned  with 
the  measurement  of  the  admittances  at  the  surfaces  of  bunting  solid  piopcllant  samples. 
The  latter  arc  needed  either  ns  inputs  m  solid  rvwket  stability  analyses  or  as  means  for 
evaluating  the  relative  driving  capabilities  of  diffeiont  solid  ptopt Hants  1 1]  While  the 
icscarch  described  herein  is  specifically  cowcctncd  with  solid  piopoll.mt  admittances  the 
developed  measurement  technique  is  quite  genet al  and  it  is  applicable  to  othci  situations 
involving  high  temperature  admittance  mcasuiements.  As  a  matter  of  fact,  the  developed 
measurement  technique  has  recently  Iven  used  successfully  in  the  measuiement  of  the 
admittances  of  gaseous  rocket  injectors  ( 2]  Other  cvamplcs  of  s<->  It  applications  include 
the  determination  of  the  impedances  of  fihrous  materials  that  ate  used  »n  exhaust  silencers 
1X4). 

Recause  of  the  complexity  of  the  combustion  process,  no  satisfactory  theoietieal 
analysis  currently  exists  for  the  accurate  determination  of  the  admittances  of  combustion 
rones  and  various  experimental  techniques  are  being  employed  to  delciminc  the  unknown 

*  hrsMtth  »>H«|s>rt«l  i««hrt  At  OXR  Grant  No  r.t.yxtl  AIAA  P,»|vr  W-tttaS.  presented  at  the  17th 
Atitivpac*  Science*  Meeting,  New  Oilcans,  USX,  Jnmmi)  IvTv 

J  Currant  tddrese  Lockheed  Georgia  Company,  tXpaiimvni  7.'-74,  Aim  *0'.  Marietta,  Ocoigta  JflOft.l, 
USA. 
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admittances.  The  development  of  such  an  experimental  technique  «»•:••  is  based  noon  a 
TsMifl  ■  >1  the  Impedance  lobe  inclh 'j  L$]  is  tl ,  '>jvi  f  .  ■ 

lu  l»M>  VtUXUVAl  impedance  tube  technique  the  experimental  arrangement  consists  01  a 
tube  with  an  acoustic  driver  at  one  end  and  the  sample  whose  admittance  is  to  be 
measured  at  the  other  end.  During  an  experiment,  the  driver  generates  an  incident  wave 
of  a  desired  frequency  that  propagates  along  the  tube  until  it  impinges  upon  the  sample 
under  test.  The  interaction  between  the  incident  wave  and  the  test  sample  results  in  a 
reflected  wave  with  modified  amplitude  and  phase.  The  reflected  wave  then  combines  with 
the  incident  wave  to  form  a  standing  wave  pattern  in  the  tube  whose  structure  depends, 
among  other  parameters,  upon  the  admittance  of  the  tested  sample.  The  structure  of  the 
standing  wave  is  measured  by  traversing  a  microphone  probe  along  the  tube.  The 
unknown  admittance  is  then  computed  by  using  the  solution  for  the  standing  wave  which 
relates  the  wave  structure  to  the  sample  admittance.  The  classical  impedance  experiment 
is  restricted  to  situations  involving  a  constant  temperature  gas  in  the  tube  and  no  mean 
llow  and  this  set-up  needs  to  be  drastically  modified  in  older  to  be  used  in  the 
mcasuicment  of  the  admittances  of  localized  combustion  piocesscs  which  introduces  a 
high  temperature  steady  flow  into  the  impedance  tube.  In  the  modified  experimental  set* 
up  (see  Figure  1)  which  was  developed  in  the  course  of  this  investigation,  the  test  sample  is 
replaced  by  the  combustion  process  under  investigation  (e  g,  the  surface  of  a  burning 
solid-propellant)  and  a  noz/lc  is  provided  at  the  ether  end  for  exhausting  the  combustion 
products.  The  acoustic  driver,  located  just  upstream  of  the  noz/lc,  is  used  to  generate  an 
incident  wave,  P,  of  a  given  frequency  /,  that  propagates  along  the  tube.  The  interaction 
of  the  incident  wave  with  the  combustion  process  icsults  in  a  reflected  wave,  PK,  with  a 
modified  amplitude  and  phase.  The  presence  of  an  axial  temperature  gradient  and  an 
axial  mean  flow  gradient  in  the  tube  icsult  in  continuous  reflection  of  the  incident  and 
reflected  waves  as  they  propagate  along  the  lube  and  the  two  systems  of  incident  ’and 
reflected  waves  interact  to  form  in  the  tube  a  standing  wave  pattern  whose  structure 
depends  upon  the  admittance  of  the  combustion  process,  the  mean  (low  properties,  the 
frequency  and  the  gas  phase  losses. 
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the  standing  wave  structure.  Instead,  the  structure  of  the  impedance  tube  standing  wave  is 
measured  by  using  several  pressure  transducers  that  arc  located  at  several  pre-selected 
locations  along  the  walls  of  the  impedance  tube.  The  measured  acoustic  pressure  data  are 
then  used  to  determine  the  wave  structure  in  the  impedance  tube  and  the  admittance  of 
the  localized  combustion  process.  In  what  follows  the  theoretical  foundation  for  the  high 
temperature  impedance  tube  will  be  described  in  detail. 


2  CONST' R VA 1  ION  KJIAIION'S 

In  this  section  the  analytical  procedure  and  numeiu.il  methods  that  have  been 
developed  to  determine  the  unknow  n  admittances  of  surfaces  exposed  to  high  temperature 
conditions  and  combustion  processes,  such  as  the  surface  of  a  burning  solid  propellant 
sample,  arc  discussed.  Specifically,  the  analytical  procedure  and  numerical  methods 
developed  to  utilize  acoustic  pressure  measurements  at  a  discrete  number  of  locations  in 
the  determination  of  the  structure  of  the  standing  wave  and  the  unknown  admittance  at 
one  end  of  the  impedance  tube,  with  account  being  taken  of  the  presence  of  an  axial 
steady  temperature  gradiQit  and  an  axial  mean  Dow  gradient  along  the  tube,  and  gas- 
phase  losses,  are  discussed. 

To  deduce  the  requited  standing  wave  struduie  from  tlu  .ncasui  emails,  one  needs  to 
obtain  solutions  of  the  conservation  equations  that  describe  the  behavior  of  small 
amplitude  distuibanccs  in  a  long  and  narrow  lube  containing  a  high  tcmpciuturc  mean 
(low  resulting  in  a  temperature  gradient  along  the  tube.  The  solutions  needed  are 
obtained  in  what  follows  here  on  the  basis  of  certain  assumptions:  (I)  the  combustion 
process  is  concentrated  at  the  surface  of  the  test  sample  in  the  tube;  (2)  the  gas  in  the  tube 
is  perfect  and  consists  of  a  single  species;  (3|  the  fluctuating  part  of  the  heat  transfer  to  the 
wall  is  negligible:  (-1)  both  the  steady  and  unsteady  flows  arc  one  dimensional:  (5)  the 
amplitudes  of  the  waves  arc  sufficiently  small  so  that  all  non-linear  leims  involving 
products  of  these  quantities  may  be  neglected:  16)  the  acoustic  gas  phase  losses  in  the  flow 
field  (e  g.,  due  to  presence  of  particles)  can  be  accounted  for  in  the  conservation  equations 
by  introducing  an  overall  bulk  loss  coefficient  or  gas  phase  loss  coefficient  which  is 
proportional  to  the  local  oscillatory  velocity.  I’ndvi  these  assumptions,  the  flow  in  the 
impedance  tube  can  be  described  by  the  following  system  of  conservation  equations  [  I  ] 


continuity: 

f  l>  <'l  i  *'(/«<)  )'  V  0. 

(I) 

momentum: 

•'(/»!<)  f'/  f  (  l/'iel  ('x  <  />  <  \  A. 

(2) 

cncigy: 

pi (i  s  i  f  i  li  i*v  <’xt 

<31 

state: 

'  p  - PRl. 

(4) 

where  p  is  the  density,  i«  the  axial  velocity,  /*  the  p'C'sme.  /  the  absolute  temperature,  v 
the  entropy,  r  the  time.  x  the  axial  co  oidinalc.  I  the  ict. tiding  foice  per  unit  xolumc 
due  to  gas  phase  losses,  Q  the  heat  gain  or  loss  pci  unit  xolumc.  R  the  specific  gas 
constant  and  C,  the  specific  heal  at  constant  xolumc 

To  obtain  the  needed  solutions,  each  of  the  dependent  variables  is  assumed  to  consist  of 
a  sum  of  a  steady  state  quantity  (to  he  denoted  by  a  Mipct posed  bar)  and  a  time 
dependent  small  amplitude  perturbation  (to  be  denoted  by  a  pnmel  These  expressions  are 
then  substituted  into  the  above  conservation  equations  which  are  then  separated  into 
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steady  state  conservation  equations  can  he  expressed  in  the  following  forms,  in  which  the 
mean  retarding  force  F  has  been  neglected: 


m*Rp  df 
m'Rf-p1  dx'* 

dp 

dx 

-U 

RT\ 

'dp  p  df  \ 
kdx  f  dx  / 

(5.6) 

du  iii  dp  _ 

mRf  / 

rdp 

pdf\ 

(7) 

dx  *  p2  dx 

P  ' 

,dx 

t  dx) 

Here  m  is  the  steady  state  mass  flow  rale  of  combustion  products  per  unit  area  that  in  the 
case  of  a  burning  solid  propellant  would  equal  p,f.  where  p,  is  the  propellant  density  and 
f  the  steady  state  propellant  burning  rate. 

Linearizing  the  corresponding  unsteady  conservation  equations,  assuming  that  the 
solutions  arc  periodic  in  time,  and  solving  for  the  space  derivatives  of  the  perturbations,’* 
one  obtains  the  following  system  of  wave  equations  [I]:  ’ 


d2, 

dx 


j 


I 

i*  i 


1.2,3, 


(8) 


where  Z,.  ll  and  Z,.  respectively,  represent  u\  p'  and  p  and  the  coefficients  AtJ  are 
A  |,  =  {up  dii/dx  -  uuup  +  Gu}/{p(r2  -  ii2)} . 


Ati  =  |(uc2/p) dp/dx  - ko}£»(c2  -  «2Jf 
Ati  *  -(I/P) dii/dx  -  {(c2ii/p) dp/dx  j£»(c2  - 
-  (i<up  4  p  dii/dx  4  C  +(M,p  dii/dx  -  it  dp/dx  t  ii»ii2p))A‘'2  -ii2), 
A}1  =  -u  {(«i2ii/p)  dp/dx  -  iruj/fc2  -  it1). 

A is  ~  u  {(c2«/p)  dp/dx}/(r2  -  a2). 


/tj,  =  -  (I,  ii|  {dp/dx  4  [«p  dii/dx  dp/dx  4  i*»r</»  t  f.u](c2  -  ii2)!, 
Ai2  »  -(l/fd{(i«-2/p) dp/dx -iru},V-*  ii2). 

.4.,  ,  =  (!,«)[  4  {(fVp)  dp/dx | /(r 2  -  ii2)]. 


(9) 

(10) 

HD 

(12) 

(13) 

(14) 

(15) 

(16) 
(17) 


where  i  =  t,p  p)* 2.  G  -  F.tt'  (e  g.,  gas  phase  loss  coefficient),  and  y  is  the  specific  heats 
ratio.  The  above  system  of  equations  represented  by  equations  |8)  can  be  considerably 
simplified  when  the  magnitude  of  the  mean  flow  u  is  very  small.  For  this  case  the 
oscillatory  flow  equations  become 

du'/dx  =  uup’lfK1.  (18) 

dp'.-dx  --  -(iujp4C)«'.  (19) 

p’  =  (p'/c2)  4  (iw'/iu) dp/dx.  (20) 


Examination  of  equations  (18)  ihrough  (.’Of  shows  that  in  this  case  the  oscillatory  density 
is  decoupled  from  the  system  of  differential  equations.  Therefore,  equations  (18)  and  (19) 
can  be  solved  (o  determine  the  oscillatory  velocity  and  pressure,  and  these  solutions 
together  with  equation  (20)  can  he  used  to  determine  the  unk now  n  density  perturbation. 

Since  the  solutions  of  the  systems  of  wave  equations  for  the  cases  6  *•  0  and 
u  *  ii|x)  *  0  are  similar,  the  method  of  solution  for  these  two  eases  can  he  illustrated  by 
discussing  (he  general  case,  when  u  P  0.  described  by  equations  (8|.  To  solve  equations  (8), 
the  coefficients  A,t  must  he  evaluated.  Inspection  of  equations  ('))  thiough  (17)  reveals  that 
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(7)  and  the  value  of  the  gas  phase  loss  coefficient  G.  Inspection  of  equations  (5)  through, (7) 
also  shows  that  the  needed  steady  state  solutions  can  be  determined  from  these  equations 
once  the  temperature  distribution  f(.\)  and  the  value  of  p(x)  at  some  x0  location  (i.e.,  the 
“initial"  point  for  the  integration)  are  known. 

Temperature  measurements  along  the  impedance  tube  could,  in  principle,  be  used  to 
provide  the  needed  axial  slcady  temperature  distribution.  However,  measured  tempera* 
ture  data  [I]  indicate  the  presence  of  a  radial  temperature  gradient  along  the  length  of  the 
tube.  Therefore,  the  question  comes  up  regarding  the  “best"  radial  location  to  be  used  for 
the  measurement  of  a  bulk  axial  temperature  distribution  that  is  consistent  with  the  one 
dimensional  wave  propagation  model  used  in  the  present  analysis.  Since  it  is  difficult  to 
locate  such  a  radial  position  in  the  impedance  tube,  a  different  approach  for  the 
determination  of  the  average  or  bulk  axial  temperature  distribution  has  been  developed. 
By  using  eneigy  conservation  considerations,  it  is  shown  in  references  [1]  and  16]  that  the 
axial  steady  state  temperature  distribution  in  the  impedance  tube  is  given  by  the  following 
expression: 

f(.v)  =  TJx)  +  (Tc  -  T.  (0))  e  *  (21 ) 

where 

A  -  C{1  +  TJx)!f[\)}\  n  -  068,  (22) 

and  the  unknown  constant  C  needs  to  be  determined  from  experimental  data.  In  the 
study  described  here  this  constant  has  been  determined  from  measured  oscillatory 
pressure  data  by  applying  a  non-linear  regression  technique,  as  will  be  discussed  below. 
The  value  of  C  thus  determined  is  the  one  that  results  in  the  “best”  fit  between  the 
measured  and  computed  impedance  tube  wave  structures.  To  complete  the  set  of  needed 
input  data,  the  mean  pressure  p(x)  can  be  measured  at  some  initial  location  x  =  x0  in  the 
impedance  tube  and  the  slcady  state  mass  flow  rate  m  cither  determined  from  known  data 
or  determined  experimentally.  For  example,  in  solid  propellant  studies  the  slcady  state 
mass  flow  rate  m  is  obtained  from  the  relation  m  =  p,r,  where  f  is  determined  from 
measured  propellant  sample  thickness  and  its  burning  time,  and  the  density  p,  is 
determined  by  weighing  a  known  volume  of  the  propellant  sample.  The  icmaining  input 
quantities,  such  as  the  flame  temperature  T(,  the  gas  constant  R  and  the  specific  heat  ratio 
y.  are  obtained  from  data  in  the  titciaiuic.  Finally,  the  wall  temperature  distribution  T. (x) 
can  be  obtained  from  measured  wall  tempo  ature  data.  With  these  data  available, 
equations  (5)  through  (7)  can  be  solved  to  dekiiuitic  the  axial  disliibution  of  all  the  steady 
slate  variables. 

Equations  (8)  are  a  system  of  three  linear,  homogeneous  coupled  firsl-oidcr  ordinary 
differential  equations,  and  they  could  be  treated  as  an  initial  value  problem  once  the 
dependent  variables  u,  p'  and  p'  were  known  at  any  x  location  in  the  impedance  tube. 
With  these  values  used  as  initial  conditions,  equations  (8)  can  be  integrated  from  the 
initial  location  toward  the  combustion  /one,  to  dctctminc  the  values  of  the  dependent 
variables  i/.  />*  and  />'  at  that  location.  Once  determined,  these  quantities  can  be  used  to 
determine  the  unknown  admittance  and  the  cntiopy  perturbation  of  the  combustion 
process. 

3.  Din  I.R  MINA  I  ION  Ob  THE  INITIAI.  CONDITIONS 

Unfortunately,  accurate  measurements  of  7,  and  /.s  (i.e.,  u'  and  (>')  arc  extremely 
difficult  if  not  impossible  at  this  time,  llcncc,  one  must  resort  to  the  use  of  a  different 
analytical  technique  for  the  solution  of  equations  (8).  Due  to  the  linearity  of  equations  (R>, 
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inS»ufedTZ<lata  in  iiie  solution  of  the  system  oi  equations  is  oased  on  Me  tiansnttstioii 
matrix"  concept  [7],  .  -  — 

Equations  (8)  can  be  written  in  the  following  matrix  form: 

{Z\  =  ML{Z»..  (23) 

where  the  subscript  x  represents  the  location  at  which  the  elements  of  the  matrix  are 
evaluated.  In  terms  of  the  backward  transmission  matrix  [?'],  {Z}  can  be  written  as 

{Z}»  -  (24) 

Equation  (24)  indicates  that  once  the  transmission  matrix  [f],  and  {Z),u  arc  known,  {Z}, 
can  be  evaluated.  To  determine  [T],  and  {Z),u.  one  can  use  equation  (24)  to  obtain  the 
follow  ing  expression  for  {Z)i(J: 

{zu  =  «5) 

Now,  differentiating  equation  (24)  with  respect  to  x  and  using  equation  (25)  yields  the* 
following  expression  for  {Z'},: 

[z\  =  [n,[n;m.  126) 

Comparing  equations  (23)  and  (26),  one  gets 

ML{Z),  =*  [TUT]X'{Z},.  (27)* 

Since  {Z}„  is  an  arbitrary  matrix,  it  can  be  eliminated  from  both  sides  of  equation  (27)  to 
give 

tn,  =  [d],[7L.  -  (28) 


The  matrix  [TJ,  can  be  determined  from  the  integration  of  equation  (28),  once  [T]tt  is 
known.  By  letting  x  go  to  x„  in  equation  (24)  one  obtains  [T],o  ={T^hc  identity  matrix. 
From  equation  (24)  it  follows  that  _ _ _  _ _ 

A.  =  in,{ZL,  Z„  =  inkAZ)„,  (29.30) 


A,  =  iri.,{Z)Mll,  Z„  =  m.  ,{*)„„•  (29.30) 

Using  equations  (29)  and  (3(T)to  obtain  expressions  for  the  elements  (Zj|M  and  (Zj),,  (eads 
to  the  derivation  of  two  algebraic  equations  involving  (Z,)>(,  and  (Zj|,„.  These  algebraic 
equations  can  then  be  solved  for  (Z,)„,  and  (2j)„,,  since  Z,  is  known  at  xn.  x,  and  x, 
from  experimental  data.  Thus,  by  making  use  of  {ZJ,,,  and  [f]„  {Zj,  can  be  evaluated 
from  equation  (24). 


4.  NON-LINKAR  Rt-GRI  SSION  TECHNIQUE 

According  to  the  above  discussion  any  three  accurate  pressure  measurements  together 
with  a  known  axial  temperature  profile  and  a  bulk  loss  coefficient  can  be  used  to 
determine  the  unknown  admittance  of  the  combustion  process.  There  arc  no  limitations 
on  the  magnitudes  of  steady  or  unsteady  variables  in  the  evaluation  of  exact  initial 
conditions  at  the  combustion  surface.  However,  errors  in  experimental  measurements 
result  in  errors  in  the  computed  admittance  values.  These  errors  can  he  minimized  by 
increasing  the  number  of  acoustic  pressure  measurements  and  utilizing  a  non-linear 
regression  technique  in  the  data  reduction  [8], 
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Let  Ei be  >  quantity  measured  »l  Xj.i  *=  1.7  .  ,«  and  7;  the  corresponding theoretically 
qu.-.-.rity  cvaiuacd  at  tl.o  same  !•  •  d  -  In  (he  pie-vnt  proNcm,  the 
experimentally  measured  quantity  is  the  acoustic  pressure  Z,  which  is  measured  at 
various  x,  locations.  Non-linear  regression  consists  of  finding  values  of  u\  />’  and  p*  at 
x  =»  x0,  (i.e.,  values  of  Z0|,  Z0J  and  Z0J)  which  give  the  best  fit  between  the  theoretically 
predicted  acoustic  pressure  distributions,  obtained  from  the  solution  of  the  differential 
equations  (23),  and  the  experimentally  measured  acoustic  pressures.  This  is  accomplished 
by  computing  the  values  of  u\  p'  and  p'  at  .x  =  x0  which  minimize  the  root-mean-square 
deviation  between  the  theoretically  predicted  acoustic  pressure  distribution  and  the 
experimental  data.  The  calculated  optimum  values  of  u\  p  and  p'  at  x  =  x0  can  then  be 
used  to  obtain  optimum  values  of  the  combustion  process  admittance  and  response 
factor.  To  find  the  minimum  root -mean-square  deviation,  the  following  function  F  is 
minimi/ed: 

F=il£|-7;i}.  (311 

i-  I 

If  a  minimum  of  F  exists,  then  the  gradient  of  F  vanishes  at  the  minimum:  that  is, 
cF/FZ0l  =  cF/cZ0,  =  rF/rZ0,  -=  0,  or 

j  =  0.  for  A  =  1.2.3.  (32) 

Equation  (32)  represents  a  set  of  three  non-linear  equations  for  the  unknowns  Z0l.  Z0l 
and  Z0j  as  T,  and  fTt/FZ0k  are  both  functions  of  ZUk .  A  Newton  Raphson  iterative 
scheme  is  used  to  obtain  a  solution  of  equation  (32),  by  utilizing  a  linearized  version  of 
this  equation  which  involves  the  expansion  of  7)  in  a  first-order  Taylor  series  with  respect 
to  the  parameters  Z0|,  Z0l  and  Z0(.  The  resulting  system  of  linear  algebraic  equations 
can  be  expressed  in  the  following  form: 

1.2,3. 
(33) 

where  m  represents  the  mth  iteration. 

For  compactness,  equation  (33)  can  now  be  rewritten  in  matrix  form,  by  making  the 
definitions 


I  {E,-T”)<'ri:Z0k 


n 

1 

u. iw  r  Z0l 

li.wlj 
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1  I 

In  terms  of  these  quantities  equation  (33)  can  be  expressed  as  3"(Z"' 1  -  H".  or. 

equivalently. 

Z,7“  =  ZJ4  M-|  ’IT.  (34) 

Equation  (31)  is  a  linear  equation  for  the  unknown  /", ' '  and  it  cun  be  readily  solved  once 
the  elements  of  Am  and  Bm  are  computed.  The  computation  of  these  elements  requires  the 
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STl  .  +a...)-ir 

r-^ok 


« 


(35) 


where  r.  is  a  small  number.  Thus,  to  solve  for  the  elements  of  the  matrices  Am  and  B"  of 
equalioh  (34),  it  is  necessary  to  compute  T,  four  times  for  each  iteration  by  solving 
equations  (23).  The  first  solution  is  obtained  at  Z~  and  the  remaining  solutions  of  T,  are 
obtained  by  varying  each  of  three  parameters  by  a  small  amount  fZ"k  in  succession. 

The  non-linear  regression  technique  can  also  be  used  to  determine  the  values  of  the 
heat  transfer  parameter  C  (i.e.,  sec  equations  (21,  22))  and  the  gas  phase  loss  coefficient  C. 
In  this  case  the  function  F  in  equation  (31)  is  partially  differentiated  with  respect  to  C  and 
C,  now  represented  by  Z0)  and  Z0j  respectively,  and  the  derivatives  arc  equated  to  zero. 
Then  all  the  steps  involved  in  the  proccduie  leading  from  equation  (32)  to  equation  (35) 
arc  repeated  to  evaluate  the  optimum  values  of  C  and  C. 

In  the  actual  computational  scheme,  the  function  F  in  equation  (311  is  expressed  cither 
hy  using  the  teal  and  imaginary  parts  of  the  oscillatoiy  pressures,  which  can  be  computed 
from  the  measured  amplitudes  and  phases,  or  by  utilizing  pressure  amplitudes  only  fl]. 
Accordingly,  the  following  functions  F  have  been  ;:scd: 


F  = 


n 


I  [</>: 

i*  i 


(36) 


or 


F  =  £<pl.,-rL>2. 

j*  I 


(37) 


where  the  subscripts  r,  i  and  a,  respectively,  represent  the  real  part,  imaginary  part  and  the 
amplitude  of  a  complex  quantity  and  the  subscript  e  refers  to  experimental  quantities.  The 
numerical  solution  of  equation  (23)  gives  the  real  and  imaginary  parts  of  the  oscillatory 
pressures  at  locations  in  the  impedance  tube  where  these  values  are  determined  exper¬ 
imentally.  Theiefore,  the  derivatives  of  these  variables  with  respect  to  u'0,  p'0  and  p'0  (i.e., 
f’p'/rZo,  and  cp',/cZ0k)  can  he  readily  computed  and  used  in  the  non-linear  regression 
scheme  where  both  the  amplitude  and  phase  of  the  oscillatory  pressures  arc  needed  to 
compute  the  admittances.  In  the  second  schcnv  where  only  the  pressure  amplitudes  are, 
used  to  compute  the  admittances,  the  derivatives  obtained  by  using 

and  f/»,V computed  in  the  first  scheme.  With  the  real  and  imaginary  parts  of  u'0.  p[, 

and  i>‘0  represented  by  ZUk,  k  —  1.2 . 6,  the  derivatives  of  pressure  amplitude  p'm  are 

given  by 


*p,  _  <\  ir,2 + p?)  _  1  ( .  <>r  .  *>;  \ 

r/.„k  "  rz„k  p.\  rz0k  PirzJ 


(38) 


Once  these  derivatives  arc  computed  the  non  linear  regression  tan  be  used  to  minimize 
the  error  between  theoretically  computed  and  experimentally  measured  pressure 
amplitudes. 


5  C  OMPUTATIONAL  St  itt  Mb  TO  M1MMI/I  I  III  I  RROR 

The  non-linear  regression  is  used  to  find  the  values  of  p\,  and  at  the  combustion 
zone  which  give  the  best  fit  between  the  ihcoictl.  ally  pr-'d icted  acoustic  pressure 
distribution  and  the  coi responding  experimentally  measured  a.  .rustic  pressures.  A  similar 
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1  MMPflBeNr  refrMSlcfh  |*I »F.iit  Is  also  "us.  !  to~  hi;iin  if ic  cons' :...W  C  .  c  these 

■"-■'tiblWB-Sncjr  rcgfCMvU  opciations  arc  coupled,  (hey  must  be  solved  together  to  obtain 
the  combustion  zone  admittance.  The  manner  in  which  this  is  done  can  be  briefly 
explained  as  follows. 

(1)  Initially,  arbitrary  values  of  C  and  G  are  chosen,  and,  with  three  measured 
pressures  used  together  with  the  “transmission  matrix”  scheme,  a  set  of  “initial” 
values  for  uj,,  p'0  and  p‘0  at  the  combustion  process  is  computed. 

(2)  The  non-linear  regression  scheme  is  then  used  to  obtain  the  values  of  Uo,  r'o  and 
and  p'0  which  minimize  the  root-mean-square  deviation  between  the  theoretically 
predicted  data  and  the  experimentally  measured  data  for  the  chosen  values  of  C 
and  G. 

(3)  Then  the  values  of  C  and  G  are  recomputed  while  keeping  (and  using)  the 
determined  (in  Step  (2))  “optimum”  values  of  Uq.  p'0,  and  pu  constant. 

(4)  Steps  (2)  and  (3)  arc  repeated  until  u’0%  p‘0 ,  p'0 ,  C  and  G  converge  to  a  given  set  of 
values. 

« 

6.  II-LUSI  RATIONS  OF  COMPUTATIONAL  SC  III  ME 

To  check  the  accuracy  of  (he  data  reduction  schemes  discussed  in  the  last  section, 
computations  were  performed  for  a  numlicr  of  hypothetical  cases  for  which  the  exact 
values  of  C.  G,  u'0,  p'0  and  p’0  wer<?  assumed  to  be  known.  With  these  values  used  as  initial 
conditions,  the  conservation  equations  (i.e.,  equations  (8))  describing  the  oscillatory  flow 
field  in  the  impedance  tube  were  solved  to  determine  the  “exact”  wave  structure  in  the 
hypothetical  impedance  tube.  The  axial  distributions  of  the  steady  state  variables  needed 
for  the  solution  of  equations  (8)  were  computed  from  equations  (5),  (6),  (7),  (21)  and  (22).  In 
these  equations  different  forms  of  wall  temperature  distributions  were  used  with  a  fixed 
flame  temperature,  Tt  of  3780°R.  Solution  of  equation  (5)  indicated  that  the  steady  state 
pressure  variation  was  almost  negligible.  Different  values  of  m  were  used  (see  equations  (5) 
_  Y  through  (7))  in  this  analysis  such  that  the  maximum  values  of  u  at  x  -  0  varied  from  0  to 
vv  200  ft/s.  The  specific  heat  ratio^v)  and  the  specific  gas  constant,  R,  were  assumed  to  be 
1-28  and  2480  ft-lbf/slug°R,  respectively. 

*  A  number  of  oscillatory  pressure  data  at  different  axial  locations  along  the  lube  were 

lakcn  from  the  solution  of  equations  (8)  and  were  perturbed  from  their  exact  values.  These 
deviations  represent  the  expected  errors  in  the  measured  acoustic  pressure  data.  These 
modified  pressure  data  were  then  Ted  into  the  developed  data  reduction  programs  in  an 
effort  to  determine  the  error  in  the  admittance  computation  resulting  from  the  artificially 
introduced  experimental  inaccuracies.  Careful  examination  of  the  results  obtained  from 
the  various  hypothetical  cases  investigated  indicates  that  the  heat  transfer  parameter  C 
converges  rapidly  to  its  correct  value,  whereas  the  gas  phase  loss  coefficient  G  converges 
slowly.  Therefore,  to  obtain  an  accurate  value  of  G  the  tolerance  limit  must  be  kept  as 
small  as  possible  and  the  errors  in  the  input  ptessure  data  should  not  be  too  large.  Two 
sets  of  such  compulations  are  presented  below,  in  these  two  cases  the  wall  temperature 
was  assumed  to  be  constant  (i.e.,  T„.  -  8KPR)  and  the  mean  flow  velocity  f»  was  kept  zero. 

Case  1.  In  this  case  the  gas  phase  loss  in  the  impedance  tube  was  assumed  to  he  zero 
and,  with  use  of  known  values  of  C,  mJj.  p'v  and  p'0,  equations  (8)  re  solved  to  determine 
the  corresponding  hypothetical  standing  wave  structure  in  the  tube.  Then,  an  arbitrary 
initial  value  of  C  was  chosen  and  input  pressure  data,  obtained  by  modifying  the 
computed  “exact"  pressure  data  at  a  discrete  number  of  points,  were  used  to  recompute 
the  wave  structure  and  the  desired  value  of  C.  The  recomputed  and  “exact"  wave 
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- structures'*  reshown  in 'Figure?  2  and  3,  and  Table  1  proxWes  *  comparison  between  the  " 

- --S#*«*.'»"“e(vV4be-com(HMc<d  value*  ,  r  ~  ,:n  /•*,  •  • 

/VlTjjure  4  shows  excellent  agioemo  >  ■  '  *  .  ;  ' 

distributions.  Therefore,  this  technique  can  be  used  to  obtain  axial  temperature  distri* 
button  in  any  pipe  with  high  temperature  flow  situations  where  direct  temperature 
measurements  are  difficult. 

Cose  2.  In  this  case  equations  (8)  were  solved  for  a  situation  where  the  gas  phase  loss 
was  assumed  to  be  present  in  the  impedance  tube  and  the  constant  G  was  assigned  a 
specific  value.  Then,  with  use  of  the  assumed  values  of  C,  G,  u‘0,  p’0  and  p'0,  equations  (8) 
were  solved  to  determine  the  hypothetical  “exact"  acoustic  pressure  wave  structure  in  the 
impedance  lube.  As  was  done  in  Case  I  above,  assumed  values  of  C  and  G  and  perturbed 
“exact"  pressure  data  were  used  to  recompute  the  wave  structure  and  the  values  of  C  and 
G.  For  brevity,  the  results  presented  herein  are  for  the  case  where/w=  01  and  they  are 
presented  in  Figures  S  through  7  and  in  Table  2.  It  should  be  pointed  out,  however,  that 
similar  results  were  also  obtained  for  wide  ranges  of  both  positive  and  negative  values^of 
the  real  part  of  the  admittance.  In  one  of  the  two  computational  schemes  used  to 
determine  the  variables  at  the  combustion  zone,  the  gas  phase  loss  coefficient  was  kept 
zero.  In  this  case  (i.e..  Plot  B  in  Figures  S  and  6)  the  recomputed  and  “exact"  wave 
structures  are  not  in  good  agreement,  particularly  at  and  near  pressure  minima,  and  there 
are  considerable  discrepancies  between  the  computed  and  exact  values  of  uj>,  p'0,  p'0  and 
the  admittance  Y.  In  the  second  case,  the  determination  of  the  bulk  loss  coefficient  G  was 
also  included  in  the  data  reduction  scheme  which  now  resulted  in  good  agreement 
between  the  two  wave  structures  (i.e.,  between  plots  A  and  C  in  Figures  S  and  6)  and  the 
variables  Uq,  p‘0.  />„  and  V.  However.  Figure  7  shows  excellent  agreement  between  the 
exact  and  both  the  recomputed  temperature  distributions  (i.e.,  between  A,  B  and  C  in 
Figure  7)  because  the  axial  steady  state  temperature  distribution  depends  only  on  the 
axial  wave  length  variation  of  the  standing  wave  and  docs  not  depend  onrthe  magnitude 
of  oscillatory  pressure. 

Since  the  primary  objective  Of  this  study  was  to  determine  whether  the  non-linear 
regression  technique  developed  could  handle  situations  in  which  a  gas  phase  loss  is 
present  in  the  flow,  the  data  reported  in  Case  2  describes  a  situation  in  which  only  small 
errors  were  introduced  into  the  “exact"  pressure  data  that  was  used  to  generate  the  input 
pressure  points  needed  for  the  data  reduction  procedure.  However,  additional,  related 
hypothetical  cases  in  which  the  input  pressure  data  contained  relatively  large  errors  were 
also  investigated  in  the  course  of  this  study  and  they  produced  good  agreements  between 
the  “exact”  and  recomputed  wave  structures  and  between  the  assumed  and  recomputed 
values  of  mJ>.  p'a.  p'0%  C ,  G  and  Y.  It  has  also  been  observed  in  these  studies  that  the  quality 
of  the  computed  results  improves  when  the  number  of  input  pressure  points  increases  and 
when  the  input  pressure  data  describe  the  wave  structure  over  a  tube  length  that  includes 
at  least  two  minimum  points. 

Examination  of  Figures  2  through  7  and  Tables  I  and  2  indicate  very  good  agreement 
between  the  assumed  and  computed  values  of  Up,  p'0,  p'0,  C,  G  and  the  admittance  V.  This 
good  agreement  suggests  that  the  data  reduction  procedures  developed  as  part  of  this 
research  effort  can  indeed  he  used  to  determine  unknown  admittance  values  in  high 
temperature  impedance  tube  experiments  without  much  loss  of  accuracy.  The  applications 
of  the  experimental  procedures  developed  in  this  paper  in  the  determination  of  the 
admittances  of  combustion  processes  are  described  elsewhere  [I.  2]. 

Reference  [1]  describes  the  application  of  this  technique  for  the  measurement  of 
burning  solid  propellant  admittances.  Some  of  the  measured  admittance  values  showed 
considerable  amount  of  scatter.  The  observed  scatter  in  the  admittance  values  was 
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Figure  3.  Exact,  recomputed  and  input  pressure  phase  data  used  to  check  the  dcxclopcd  data  reduction 
procedure.  )'  «  - 005.  C  •>  0.  chamber  pressure  -  15  7  psi,  excitation  frequency  --  HKIO  H<;  .  exact 

structure;  -  -  .  structure  recomputed  by  using  pressure  amplitudes  and  phases:  .  structure  recomputed 

by  using  pressure  amplitudes  only :  •.  "perturbed'  input  pressure 
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Figure  4.  Axial  distribution  of  exact  and  recomputed  steady  state  temperature  in  the  impedance  tube 
xc-tRB  C*»*0  1  — kai  pr.neiini  f  7  mi.  cxi-.f ,| I i, >11  fii  n III  nn  — l‘ni|  .  exact  structure; 

structure  recomputed  by  using  pressure  amplitudes  and  phases.  -  .  stiucturc  iccomputed  by  using  pressure 

amplitudes  rnlj ;  #i  Ttrtrrhcil"  input  pr^sum-  smpi . .  jhre 
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Comparison  IvlHvcfi  cc.wf  ttdues  and  <  •  %  ■•••  ,• 

technique 


Chamber  pressure  =  15  7  psi 
Frequency  of  excitation  -  1000  Hz 


• 

Parameters 

Collect 

values 

Computed 
using 
amplitude 
and  phase 

Computed 

using 

amplitude 

only 

Initial 

guess 

Rea)  part  of  admittance.  V, 

-005 

0  049757 

0052565 

Imaginary  part  of  admittance.  )’ 

015 

01236 

01232 

Heat  transfer  paiumcier,  C  (ft) 

0943 

0972 

0973 

2-0 

Pressure  at  Amplitude  (d B) 

145436 

145815 

145  SIS 

* 

combustion  zone  p'a  Phase  (deg) 

2302 

2299 

2286 

Velocity  at  Amplitude  (ft/s) 

148 

1  301 

1  308 

* 

combustion  zone  Phase  (deg) 

-214 

-1817 

-18  27 

Tensity  at  Amplitude  (slug/l’)  0612  x  10*  0-6472x10  *  0  6477x10'* 

combustion  zone  p'0  Phase  (deg)  /  229  0  228-7  2274 


Ta«i  r  2 

Data  computed  with  g as  phase  loss ,  by  using  non-linear  regression  for  a  positive  real 

admittance  * 

Chamber  pressure  -  300  psi 
Frequency  of  excitation  -1000  Hz 


Computed  values  using 
amplitude  and  phase 


Loss  coefficient 

and 

heat  transfer 

parameter 

Correct 

Assuming 

computed 

Pm  a  meters 

values 

zero  loss 

simultaneously 

Real  part  of  admittance,  V, 

01 

0055532 

0)0147 

Imaginary  part  of  admittance.  >', 

015 

013315 

01362 

Rulk  loss  coefficient.  O'  (fbf  s,ft*| 

15 

OO 

1556 

Heat  transfer  parameter,  C  (ft) 

10 

1024 

1021 

Pressure  at 

Amplitude  (dB) 

1454357 

1454 

1455 

combustion  zone 

Phase  (deg) 

2301944 

2299 

2301 

Velocity  at 

Amplitude  (ft/s) 

0088204 

007068 

008356 

combustion  zone  Wq 

Phase  (deg) 

75496 

62  76 

-  76  54 

Density  at 

Amplitude  (slug  ti’) 

06143 x  10  * 

06197x10  T 

06219x10  * 

combustion  zone  p‘„ 

Phase  (deg) 

2324994 

231  1 

2324 
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figure  5  F xact,  recomputed  and  inpul  pressure  amplitude  data  used  in  check  tin-  developed  data  reduction 
procedure;)',  -  0 1.  G  -  1  5.  chamber  pressure  -  W)psi,  excitation  frequency  lIXtOH/;  .  exact  structure 
(A);  recomputed  structure  when  overall  hulk  loss  coefficient  was  V ^ pi  zero  (B).  -  recomputed 

structure  when  overall  bulk  loss  was  considered  to  he  present  (C'l;  #.  "1*11111  bed"  input  pressure  amplitude 
data.  ’ 


figure  6.  F*act,  recomputed  and  input  pressure  phase  data  used  to  check  the  developed  data  reduction 
procedure;  Y,  =  0  I,  G  =  I  S.  chamber  pressure  =  300  psi,  excitation  frequency  --  1000  Hz;  .  react  sli  uclurc 
(Al.  -  •  recomputed  structure  when  cserall  bulk  loss  coefficient  was  kept  zero  (Hi;  ,  recomputed 

structure  when  overall  hulk  loss  was  considered  lo  be  present  (C|;  9.  "perturhed"'  input  pressure  phase  data 


Figure  7.  Aval  dislribui  r  n  ofesaei  and  recomputed  steady  stale  rcmpcr.nure  in  the  miprdanc c  lube; 
exact  structure  (Al;  .  recomputed  structure  when  overall  bulk  loss  coefficient  was  kept  zeto  (Bl; 
recomputed  structure  when  mctall  bulk  loss  was  considered  to  be  present  (C'l. 
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describe  the  standing  wave  in  the  impedance  tube.  However,  in  a  parallel  study  [2]  where 
these  shortcomings  were  overcome  with  a  steady  burning,  longer  test  duration  and  large 
number  of  pressure  measurement  points  covering  two  or  more  pressure  minima  to 
describe  the  standing  wave,  the  admittance  values  obtained  showed  very  little  scatter  and 
they  agree  very  well  w  ith  the  theory.  In  this  study  [2]  the  admittances  of  reactive  gaseous 
rocket  injectors  were  measured  by  using  the  modified  impedance  tube  technique. 
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INTRODUCTION 


This  paper  describes  recent  advances  made  in  a  research  program,  sponsored 
by  AFOSR,  that  is  concerned  with  the  application  of  the  modified  impedance  tube 
technique  in  the  measurement  of  the  admittances  and  response  functione  of 
burning  solid  propellents.  This  data  is  needed  to  evaluate  the  relative  driving 
characteristics  of  various  solid  propellants  and  as  an  input  data  for  combustion 
instability  analysis  of  solid  rocket  motors.  To  date,  most  response  function 
measurements  were  performed  utilizing  the  T-burner  technique'  ' •  However,  due 
to  difficulties  associated  with  the  interpretation  of  T-burner  data,  additional 
experyng^tal  techniques,  such  as  the  rotating  valve'  '  and  the  impedance 
tube'  ’  ,  have  recently  been  used  to  measure  propellant  response  functions. 

A  schematic  of  the  modified  impedance  tube  is  shown  in  Fig.  1.  The  tested 
solid  propellant  disc  is  placed  at  one  end  of  the  tube  and  a  combination  of  an 
acoustic  driver  and  an  exhaust  valve  is  placed  at  the  opposite  end.  Performing 
a  test  consists  of  turning  on  the  acoustic  driver  to  excite  a  standing  wave  of  a 
predetermined  frequency  in  the  tube  and  the  ignition  and  burn  out  of  the  teated 
solid  propellant  sample.  During  the  test,  heat  losses  to  the  tube  walls  and  the 
presence  of  waves  moving  in  both  directions  result  in  an  axially  varying  steady 
temperature  profile  and  a  standing  acoustic  wave  in  the  tube,  as  shown  in  ««. 
2,  It  can  be  shown  chat  the  characteristics  of  the  standing  wave  are  dependent, 
among  other  things,  upon  the  boundary  conditions  at  the  burning  propellant 
surface  that  are  described  by  specifying  the  admittance  at  the  propellant 
surface.  It  can  be  shown v  that  the  magnitude  and  sign  of  the  real  part  of  the 
admittance  are  directly  relatod  to  the  magnitude  and  direction  of  acoustic 
energy  flow  at  the  combustion  zone.  When  the  real  part  of  the  admittance  is 
positive,  acoustic  energy  is  flowing  from  the  combustion  sons  into  the  tube  and 
the  propellant  is  said  to  be  driving.  When  this  occurs,  disturbances  in  the  gas 
phase  are  amplified  by  the  combustion  process  at  the  propellant  surface. 

In  the  impedance  tube  experiment,  the  relationship  between  the  admittance 
of  the  burning  solid  propellant  surface  and  the  impedance  tube  wave  structure  is 
'•r.iised  to  determine  the  admittance  at  the  propellant  surface.  This  is 
accomplished  by  measuring  the  resulting  acoustic  wave  structure  and  relating  it 


*  Research  Supported  under  AFOSR  Grant  73-2571 

*•  Graduate  Research  Assistant 

#  teeier  Research  Engineer 

M  Regents*  Professor 


I-A-18 


to  th«  unknown  admittance  by  utilizing  the  solutions  of  the  impedance  tube  wave 
equations. 

Fifteen  pressure  transducers  are  placed  at  predetermined  locations  along 
the  tube  walls  to  measure  the  acoustic  pressure  amplitudes  and  phases.  The  data 
acquired  at  different  instances  during  the  test  period  are  input  into  a  newly 
developed  data  reduction  scheme  that  determines  the  admittances  and  response 
functions  at  the  burning  solid  propellant  surface.  This  paper  describes  recent 
efforts  that  were  concerned  with:  (1)  the  improvement  of  the  data  acquisition 
system;  and  (2)  the  development  of  a  new  data  reduction  scheme. 

IMPROVED  DATA  ACQUISITION  SYSTEM 

During  the  past  several  months,  efforts  have  been  made  to  improve  the 
accuracy  of  the  measured  data.  The  efforts  included:  (1)  increasing  the  accuracy 
and  the  dynamic  range  of  the  measured  acoustic  data  and  (2)  increasing  the 
number  of  pressure  transducers  utilized  in  the  determination  of  the  acoustic 
wave  structure  in  the  impedance  tube.  A  combination  of  sixteen  piezoelectric  and 
condenser  type  pressure  transducers  wore  utilized.  To  acquire  the  pressure 
amplitude  and  phase  data  measured  by  the  ten  piezoelectric  pressure  transducers, 
a  minicomputer  data  acquisition  system  is  utilized.  This  system  consists  of  a 
Hewlett  Packard  2100S  minicomputer  with  an  HP  7901  disc  system  and  a  Preston 
GMAD-1  analog-to-digital  convector.  By  passing  the  signals  measured  by  the 
pressure  transducers  through  low  pass  filters,  the  D.C.  component  is  filtered 
out  and  the  A.C.  component  is  amplified  to  levels  below  the  A-to-D  input  voltage 
limitation.  The  software  utilized  for  data  acquisition  via  this  system  has  been 
rewritten  to  include  direct  memory  access,  resulting  in  an  increase  of  the 
maximum  transfer  rate  to  the  minicomputer  from  40,000  words  per  second  to 
600,000  words  per  second. 

In  addition,  theoretical  studies  of  the  characteristics  of  the  standing 
wave  structure  have  indicated  that  the  accuracy  of  the  measured  admittances 
depends  upon  the  ability  to  accurately  measure  the  difference  between  the  maxima 
and  minima  of  the  standing  wave  amplitudes;  To  acquire  such  a  capability,  the 
transducer  calibration  procedure  has  been  modified.  Transducers  located  near 
the  pressure  minima  are  calibrated  to  cover  the  range  110-155  dB  while 
transducers  located  near  pressure  maxima  arc  calibrated  to  cover  the  120-165  dB 
range,  as  shown  in  Fig.  3.  This  procedural  results  in  a  measurement  system 
capable  of  measuring  a  110-165  dB  range  that  is  needed  to  obtain  the  desired 
experimental  accuracy. 

Previous  analytical  studies  had  indicated  that  reliable  determination  of 
the  admittance  of  a  burning  solid  propellant  require  the  determination  of  the 
standing  wave  structure  through  a  distance  covering  at  least  two  standing  wave 
minima.  This  is  illustrated  in  Fig.  4  which  describes  a  hypothetical, 
theoretically  predicted  standing  wave  in  an  impedance  tube  for  a  given 
admittance  value  and  different  values  of  the  gas  phase  loss  parameter  G  that  is 
representative  of  the  acoustic  energy  loss  in  the  gas  phase.  Fig.  4  shows  that 
pressure  amplitudes  corresponding  to  two  different  values  of  G  coincide  along  a 
distance  measured  from  the  propellant  sample  (i.e.,  x  •  0)  through  the  first 
pressure  minimum,  but  they  "separate  out"  as  the  second  minimum  point  is 
approached.  This  clearly  indicates  that  for  an  accurate  measurement  of  the  bulk 
loss  parameter  G,  and  hence  the  admittance  of  the  tested  solid  propellant,  the 
standing  wave  structure  has  to  be  measured  through  a  distance  covering  at  least 
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two  standing  wave  pressure  minima.  To  accomplish  this,  the  mnber  of  pressure 
transducers  has  been  increased  from  ten  to  fifteen.  A  typical  distribution  of 
the  pressure  transducers  along  the  tube  is  shown  in  Pig.  5.  It  is  worth  noting 
that  the  more  accurate  piezoelectric  pressure  transducers  are  located  near 
pressure  minima  where  the  greatest  experimental  accuracy  is  required. 

IMPROVED  DATA  REDUCTION  PROCEDURE 


Earlier  in  this  program,  the  impedance  tube  wave  equations,  shown  in  Fig. 
6,  were  utilized  in  an  iterative  solution  procedure  to  determine  in  the  unknown 
admittance  at  the  propellant  surface  and  the  parameters  G  and  C  that 
respectively  described  the  gas  phase  losses  and  the  steady  state  temperature 
distribution.  Subsequently,  this  procedure' was  improved  by  eliminating  the  need 
to  numerically  iterate  while  determining  the  unknown  propellant  admittance.  The 
latter  change  resulted  in  considerable  savings  in  computer  time. 

More  recently,  the  accuracy  of  the  determined  unknowns  haa  been  improved 
by  redefining  the  error  function  F  whose  minimization  results  in  the 
determination  of  the  needed  unknowns.  This  improvement  is  outlined  in  Fig.  7 
where  both  old  and  newly  defined  error  functions  are  shown.  When  the  previous 
error  function  was  utilized,  errors  near  pressure  amplitude  maxima  completely 
dominated  the  behavior  of  the  error  function  with  errors  near  minisia  point 
having  little  influence  upon  the  behavior  of  F.  It  shows  in  Fig.  7  that  when  the 
previous  definition  of  F  was  used,  a  1  dB  error  near  a  maximum  point  resulted  in 
a  "contribution"  to  F  that  is  1000  times  larger  than  a  corresponding 
contribution  from  a  1  db  error  near  a  minimum  point.  Also  shown  in  Fig.  7  is 
the  newly  defined  error  function  F  where  a  "weighting"  procedure  is  used  to 
eliminate  the  discussed  uneven  contributions  to  F  of  different  measurement 
locations. 


RESULTS 


Considering  the  complexity  of  the  propellant  admittance  measurement 
technique,  the  repeatability  of  the  experimental  data  was  of  snich  concern.  The 
repeatability  achieved  is  demonstrated  in  Fig.  8  showing  the  time  evolution  of 
pressure  amplitudes  measured  at  distances  of  2.0  and  11.75  inches  from  the 
propellant  surface  in  two  different  tests  conducted  with  an  A-13  propellant. 
Examination  of  this  figure  indicates  remarkable  similarity  in  the  time  variation 
of  the  amplitudes  and  excellent  agreement  between  the  computed  values  of  the 
admittance  and  the  bulk  gas  phase  loss  G. 

Typical  admittance  and  gas  phase  loss  data  obtained  during  tests  conducted 
with  the  composite  propellants  A-13  and  A-14,  at  300  psig  chamber  pressure  are 
presented  in  Figures  9  through  13.  Figures  9  and  10  show  the  time  variations  of 
the  real  and  imaginary  parts  of  the  admittance  Y  and  the  bulk  gas  phase  loss  G 
during  tests  conducted  with  A-13  and  A-14  propellants.  It  is  noted  that  the  real 
part  of  the  admittance  YR  varies  little  with  time  while  energy  losses  in  the 
gas  phase  vary  significantly  with  time.  These  figures  demonstrate  the  ability  of 
the  impedance  tube  technique  to  simultaneously  measure  the  time  variation  of 
acoustic  energy  gains  and  losses  inside  the  tube  during  a  test.  Figures  11  and 

12  describe  the  frequency  dependence  of  the  real  part  of  the  admittance  Yo  and 
the  bulk  gas  phase  loss  parameter  G  obtained  from  experiments  conducted  with  A- 

13  and  A-14  propellants.  The  YR  curves  indicate  that  this  quantity  peaks  at 
around  650  Hz  for  A-13  propellant  and  at  around  1000  Hs  for  the  A-14  propellant. 
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A  comparison  between  values  of  the  real  part  of  the  response 
function  R^'1^  »  for  ,the  A-13  propellant,  obtained  ,by  the  impedance  tube 
technique,  the  T-burnerU  ,  and  the  microwave  technique'07  is  shown  in  Pig.  13. 
Good  agreement  is  demonstrated  at  frequencies  below  450  Ha  while  at  frequencies 
above  450  Ha  the  results  differ  significantly.  The  reasons  for  these 
discrepancies  are  yet  to  be  resolved. 

CONCLUSIONS 


t 


Significant  improvements  have  been  achieved  in  the  data  acquisition  and 
data  reduction  procedures,  as  demonstrated  by  the  results  provided  in  this 
paper.  It  has  been  shown  that  the  real  part  of  the  admittance  remains 
constant  during  the  steady  state  burn  period  while  gas  phase  losses  may  vary 
significantly  during  the  same  time  period.  Reproducibility  of  the  measured 
quantities  has  been  established  as  well  as  the  capability  of  this  experimental 
technique  to  measure  gas  phase  losses  in  the  impedance  tube  simultaneously  with 
the  measurement  of  the  burning  solid  propellant  surface  admittance. 
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An  Impedance  Tube,*  Typical  Axial  Steady  State  Temperature 
Distribution  and  a  typical  Standing  Wave  Pattern. 
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In  obtaining  a  solution  the  above  equations  are  recast  in  the  following 

form 
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Where  f  j.  fj  and  f ^  are  functions  of  u,  p/,  p',  c ,  C  and  «. 
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Fig.  6.  Impedance  Tube  Wave  Equations. 
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Modified  Data  Reduction  Scheme. 
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Fig*  11*  Real  Part  of  the  Admittance  vs.  Frequency  A-13  and  A-14 
Propellant*  300  p*ig. 
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Fig.  12.  BULK  LOSS  COEFFICIENT  VS.  FREQUENCY 
A- 1 3  and  A- 1 4  Propellants,  300  psig 
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DETERMINATION  OP  SOLID  PROPELLANT  ADMITTANCES  BY  THE 
IMPEDANCE  TUBE  METHOD* 
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School  of  Aoroapoco  Engineering 
Georgia  Institute  of  Technology 
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Abatract 

This  paper  describes  the  adaptation  el  the 
impedance  tube  technique  for  the  measurement  of  solid 
propellant  admittances  and  response  functions.  These 
quantities  are  needed  for  combustion  stability  analysis  of 
solid  rocket  motors.  The  experimental  set  up  consists  of  a 
tube  with  a  disc  of  solid  propellant  sample  placed  at  one 
end  and  a  combination  of  an  exhaust  valve  and  an  acoustic 
driver  placed  at  the  other  end.  Performing  a  test  consists 
of  turning  on  the  acoustic  driver  to  excite  a  starring  wave 
of  a  predetermined  frequency  in  the  tube,  ignition  and 
burnout  of  the  solid  propellant  sample.  During  a  test, 
acoustic  pressure  data  is  measured  by  fifteen  pressure 
transducers  distributed  along  a  distance  that  includes,  at 
least,  two  standing  wave  pressure  minima.  The  data  is 
transferred  via  a  fast  analog-to- digital  convertor  to  a  mini 
computer  system  for  immediate  storage  and  later  analysis. 
The  measured  acoustic  pressure  amplitude  and  phase  data 
are  input  Into  a  newly  developed  data  reduction  procedure 
that  is  based  upon  the  solution  of  the  impedance  tube  wave 
equations,  for  the  determination  of  the  admittance  and 
response  function  of  the  tested  solid  propellant  sample  and 
the  acoustic  energy  losses  in  the  gas  phase.  The  paper 
presents  results  obtained  in  a  series  of  tests  conducted  with 
an  A- IS  propellant.  The  paper  demonstrates  the  capability 
of  the  developed  experimental  technique  to  simultaneously 
determine  the  acoustic  characteristics  of  the  flow  inside 
the  impedance  tube  and  the  admittance  of  the  burning  solid 
propellant  for  the  duration  of  the  experiment. 

Nomenclature 

A-  coefficients  defined  in  Ref.  7 

C  a  constant  defined  in  Eq.  < 

c  velocity  of  sound,  in/sec. 

Cp  specific  heat  at  constant  pressure,  BTU/Lbm,  R 

cy  specific  heat  at  constant  volume,  BTU/Lbm, °R 

E  error  function  defined  in  Eq.  IS  (Lbf/in2)2 

F  losses  due  to  viscosity  and  gas  phase  damping, 

Lbm/insec 

G  gas  phase  bulk  loss  coefficient 

I  identity  matrix 

M  Mach  number,  dimensionless 

m  mass  flaw  rate  per  unit  cross  sectional  area, 

Lbm/in  - 

p  pressure,  Lbf/ln* 

0  volumetric  heat  source 

R  specific  gas  constant  Lbf-ft/slug  R  ,  also  solid 

propellant  response  function,  timenslonless 
r  burning  rate,  in/sec 

s  entropy,  BTU/LbnrR 

T  temperature,  Tt  i  also  Transmission  matrix 

defined  in  Eq.  10 

T  wall  temperature  at  x  ■  0,  R 
o 
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u  axial  velocity,  in/sec 

Y  specific  admittance,  dimensionless 

Z  variable  used  to  represent  oscillatory  quantities 

(defined  in  Eq.  6)  , 

p  density  of  solid  propellant,  Lbm/in 

a*  density,  Lbm/in* 

Superscripts 

n  variable  describing  a  steady  state  quantity 
(  )'  variable  describing  a  perturbation  quantity 

Subscripts 

(  >b  quantity  evaluated  at  the  propellant  surface 
(  )  quantity  evaluated  at  the  wall 

(  y  real  part  of  a  complex  quantity 
(  ^  imaginary  part  of  a  complex  quantity 


Introduction 

This  paper  describes  recent  advances  made  in  the 
application  of  the  modified  Impedance  tube  technique  in  the 
measurement  of  solid  propellant  admittances  and  response 
functions.  The  solid  propellant  admittance,  defined  as  the 
complex  ratio  of  the  velocity  perturbation  normal  to  the 
propellant  surface  and  the  pressure  perturbation  evaluated 
at  the  proepltant  surface.  Is  a  measure  of  the  amplification 
(or  attenuation)  sustained  by  a  gas  phase  disturbance  upon 
interaction  with  the  combustion  process  at  the  solid 
propellant  surface.  Consequently,  the  admittances  of  solid 
propellants  are  used  to  evaluate  the  relative  driving 
characteristics  of  various  solid  propellants.  In  addition,  since 
the  stability  of  a  solid  rocket  depends  upon  the  energy 
balance  between  various  processes  that  contribute  to  the 
combustor  disturbance  attenuation  or  amplification,  the 
ability  to  determine  quantitatively  the  admittance  of  the 
burning  solid  propellant  it  of  utmost  importance  for  the 
development  of  a  capability  for  designing  stable  rocket 
motors.  The  development  of  an  experimental  technique 
capable  of  measuring  the  admittance  and  response  functions 
of  burning  solid  propellants  under  conditions  simulating 
those  experienced  by  unstable  solid  rockets  is  the  main 
objective  of  the  research  work  described  in  this  paper. 

To  date,  most  response  function  measurements  were 
performed  utilizing  the  self  excited  T-bumer  technique'  \ 
When  weakly  driving  and  aluminized  propellants  failed  to 
generate  the  desired  .oscillations  in  the  T-burner,  the 
variable  area  T-bumer1  and  the  pulsed  T-bumer  were 
developed.  However,  due  to  difficulties  associated  with  the 
Interpretation  of  T-bumer  data,  adftttonal  experimental 
techniques  such  as  tb^  rotating  valve  ,  microwave  and 
the  impedance  tube'®’  ’  have  recently  been  developed  for 
the  measurement  of  burning  solid  propellant  response 
functions.  Comparisons  of  the  data  generated  by  these 
various  experimental  techniques  would  hopefully  result  tat 
accurate  solid  propellant  response  function  data. 

A  schematic  of  the  modified  impedance  tube  is  shewn 
tat  Fig.  1.  The  tested  solid  propellant  dUc  is  placed  at  one 
end  of  the  tube  and  a  combination  of  an  acoustic  driver  and 
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an  exhaust  valve  it  placed  •»  At  opposite end. perform¬ 
ing  a  test  consists  of  turnip  on  the  acoustic  driver  to 
tacit*  a  stendtog  ««  of  a  preAtarmtood  frequency  to 
At  tube  and  At  ignition  and  bum  out  of  A*  totted  solid 
propellent  sample.  Ptrtg  Ae  test,  heat  totttt  to  At  btot 
watt*  Md  At  presence  of  waves  mottos  to  both  (traction* 
result  to  to  axially  varying  tttody  temperature  profUt  and 
a  atendtos  acoustic  wave  in  At  tub*.  It  can  bt  shown  that 
At  characteristics  of  At  resulting  standtog  wave  are 
Apendmt,  among  other  things,  upon  the  boundary 
cwdltions  at  the  burning  propellant  surface  that  can  be 
described  by  specif wtog  the  propellant  surface  admittance. 
K  can  also  be  shownwrthat  the  magnitude  and  sign  of  the 
real  part  of  the  admittance  are  durectiy  related  to  Ae 
magnitude  and  direction  of  acoustic  energy  flow  at  the 
bumtog  propellant  surface.  When  the  real  part  of  the 
admittance  is  positive,  acoustic  energy  is  flowing  from  Ae 
combustion  zone  into  the  impedance  tube  and  the 
propellant  is  said  to  be  driving.  When  this  occurs, 
disturbances  in  the  gas  phase  are  amplified  by  the 
combustion  process  at  Ae  propellant  surface. 
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In  the  Impedance  tube  experiment,  the  relationship 
between  the  admittance  at  the  burning  solid  propellant 
surface  and  the  impedance  tube  wave  structure  is  utilized 
to  determine  the  admittance  at  the  propellant  surface. 
This  is  accomplished  by  measuring  At  resulting  acoustic 
wave  structure  and  relating  it  to  Ae  uiknown  admittance 
by  utilizing  the  solutions  of  the  impedance  tube  wave 
equations. 

During  a  test,  fifteen  pressure  tranedueers  are 
placed  at  predetermined  locations  along  the  tube  trails  to 
measure  Ae  acoustic  pressure  amplitudes  and  phases.  The 
Ate  acquired  at  Afferent  instances  during  the  test  period 
are  Input  into  a  newly  developed  data  reduction  scheme 
that  determines  the  admittance  and  response  function  at 
Ae  burning  solid  propellant  surface.  This  piper  describes 
recent  efforts  thot  were  concerned  wlA  improvement  of 
the  Ate  acquisition  system  and  Ae  development  of  a  now 
dan  reduction  scheme,  to  addition,  results  obtained 

•  Under  isenfropic  conditions,  otherwise  Ae  entropy 

perturbation  of  the  surface  to  also  needed 
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to  tests  conducted  wiA  a  composite  solid  propellant  (Le* 
A- IS)  are  presented  and  dtocuseed. 

Improved  Date  Reduction  Procedure 

A  new  Ate  reduction  procedure  that  utilises 
acoustic  pressure  amplitude  and  phase  Ata  measured  at 
dtocrete  points  along  the  wall  of  the  impedance  tube  for 
the  determination  of  the  admittance  of  the  tented  solid 
propellant,  to  outlined  in  this  section.  To  Avolop  the 
needed  Ate  reduction  scheme,  the  assumptions  introduced 
to  Ref.  (7)  for  describing  the  impedance  tube  flow  are 
utilized.  Under  these  assumptions,  the  flow  inside  the  tube 
can  be  described  by  the  following  set  of  one  dimensional 
conservation  equate  * 


Continuity 

£♦&<*«>  -o 

(1) 

— - 

£<P«)+£<puJ>  -*£-f 

(2) 

Energy 

**<  H  +  *  A  "  si 

(3) 

State 

P  »  p*T 

<*> 

First  and  second  laws  of  thermodynamics! 
ds  -  cp  S 

(5) 

Next,  each  of  the  following  dependent  variables  is 
expressed  as  a  sum  of  a  tteady  state  tolution  and  a  time 
dependent  perturbation  (e.g.,p(x,t)»p(x)4f>'(x,tftaiidAese 
expressions  are  substituted  into  Eqs.  (1)  through  (3). 
Separating  the  resulting  equations  into  the  corresponding 
steady  state  and  unsteady  systems  of  equations,  linearizing 
the  unsteady  equations  assuming  that  the  solutions  are 
periodic  in  time  (i.e.,p'(x,t)aelwt),  as  is  the  case  in 
the  impedance  tube  experiment),  using  the  relation  F'  * 
Gu'  (where  G  is  a  constant)  yield  a  system  of  linear  wave 
equations  that  can  be  expressed  in  the  following  form 


Where  Zj  ,  Z2  and  2  3  respectively  represent  u  ,  p  and  o'. 
The  coefficients  A  ij  depend  upon  the  corresponding 
steady  state  solutions  and  theacoustic  energy  gas  phase  loss 
that  to  specified  by  e  parameter  G.  Expression  for  Ae 
coefficients  Ajj  are  given  in  Ref.  7. 

Before  proceeding  with  the  solution  of  Eqs.  (6),  a 
comment  regarding  the  determination  of  the  steady  state 
solution  to  in  order.  Analysis  of  the  correapondtog  steady 
state  conservation  equations  shows  that  once  the  behavior 
of  one  of  the  dependent  variables  to  known,  the  remaining 
dependent  variables  could  A  determined  by  utilizing  the 
steady  state  conservation,  equations.  Experiments  and 
related  analytical  studies1'  conducted  earlier  under  this 
investigation  had  shown  thot  a  steady  state  temperature 
dtotribution  that  to  consistent  wiA  the  one  dimensional 
Ascription  of  the  wave  propagation  problem  to  the 
impedance  tube  to  provided  by  Ae  following  expression! 


*<*>-y*>  +  (Te-  Tf(o»o‘*/A 
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yields  a  unfoye  set  of  initial  cendtitone  u,fl  andp  at  the 
buniini  propellant  surface  (U.(  *  «  0)  that  minimizes  tha 
error  (unction  E. 


Determlnation  at  tha  Prooeilant  Admittance 
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calibrated  to  covor  tha  1 20- 14  J  dB  range,  at  dtown  in  Pig. 
2.  Thia  procedural  reaulu  in  a  measurement  system  capable 
of  measuring  a  1 10-149  dB  range  that  la  needed  to  obtain 
tha  daalrad  experimental  accuracy. 


Tha  following  procedure  has  been  utUisad  in  the 
determination  of  the  tmknown  propellant  admittance  from 
the  measured  acoustic  preseive  da  tat 

(1)  Initially,  arbitrary  value*  of  G  and  C  are  assumed. 
Next,  [A  J*  is  evaluated.  Equations  (12)  are  then 
bttegrated  forward  utilizing  a  fourth-order  Runga-Kutta 
method  to  evaluate  [Tj* 

(2)  Utilizing  die  calculated  [T]x  and  the  set  of 

acoustic  pressure  measurements  P*.  ,  the  values  of  u^p' 

and#  are  evaluated  at  the  propellant  surface  by  utilizing 
Eq.  (21 U 

(3)  [ay*  calculated,  utilizing  Eq.  (10). 

(t)  E  Is  calculated,  utilizing  Eq.  (19). 


(9)  The  above  calculated  error  function  E  represents  its 
minimum  for  a  given  set  of  the  parameters  C,  G  diet 
appear  in  the  matrix  L  A  j*,To  minimize  the  error  function 
E  with  respect  to  G  and  C,  Muller's  iteration  technique  has 
been  utilized.  This  technique  provides  a  rapid  method  for 
alternately  iterating  on  the  parameters  G  and  C  until  the 
minimum  of  E,  lor  which  ]§•  ■  if-  *  0  »  ,  is 

reached.*  The  admittance  and  the  response  function  of  the 
burning  solid  propellant,  the  g as  phase  acoustic  energy 
losses  and  the  axial  distribution  of  the  acoustic  velocity, 
density  and  entropy  are  evaluated  for  the  minimizing  set  of 
fr^GandC. 

Improved  Data  Acquisition  Procedure 

During*  the  past  several  months,  efforts  have  been 
made  to  improve  the  accuracy  of  the  measured  data.  The 
efforts  include*  (1)  increasing  the  accuracy  and  the 
dynamic  range  of  the  measured  acoustic  data  and  (2) 
increasing  the  number  of  pressure  transducers  utilized  in 
the  determination  of  the  acoustic  wave  structure  in  the 
impedance  tube.  A  combination  of  fifteen  piezoelectric 
and  condenser  type  pressure  transducers  were  utilized.  To 
acquire  the  pressure  amplitude  and  phase  data  measured  by 
the  ten  piezoelectric  pressure  transducers,  a  minicomputer 
data  acquisition  system  is  utilized.  This  system  consists  of 
a  Hewlett  Packard  2100S  minicomputer  with  an  HP  7901 
tfisc  system  and  a  Preston  CM  AD- 1  analog-to-dlgital 
convertor.  By  passing  the  signals  measured  by  the  pressure 
transducers  through  high  pass  filters,  the  d.c.  component  is 
filtered  out  and  the  a^.  component  is  amplified  to  levels 
below  the  A-to-D  input  voltage  limitation.  The  software 
utilized  for  data  acquisition  via  this  system  has  been 
rewritten  to  include  (firect  memory  access,  resulting  in  an 
increase  of  the  maximum  transfer  rate  to  the 
minicomputer  from  SO, 000  words  per  second  to  400,000 
words  per  second. 

In  addition,  theoretical  studes  of  the  characteristics 
of  the  standing  wave  structure  have  indkated  that  the 
accuracy  of  the  measured  admittances  depends  upon  the 
ability  to  accurately  measure  the  difference  between  the 
maxima  and  minima  of  the  standtag  wave  amplitudes.  To 
acquire  such  a  capability,  the  transducer  calibration 
procedure  has  been  modified.  Tranducers  located  near  the 
pressure  minima  are  calibrated  to  cover  the  range  110-199 
dB  while  transducers  located  near  pressure  maxima  are 
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Previous  analytical  studies  had  indicated  that 
reliable  determination  of  the  admittance  of  a  burning  solid 
propellant  requires  the  determination  of  the  standing  wave 
structure  through  a  distance  covering  at  least  two  standing 
wave  minima.  This  is  illustrated  in  Pig.  3  which  describes  a 
hypothetical,  theoretically  predicted,  standing  wave  in  an 
impedance  tube  for  a  given  admittance  value  and  different 
values  of  the  gas  phase  loss  parameter  G  that  is 
representative  of  the  acoustic  energy  tosses  in  the  gas 
phase.  Figure  3  shows  that  pressure  amplitudes 
corresponding  to  two  different  values  of  G  coincide  along  a 
distance  measured  from  the  propellant  sample  (i.e.,  x  «  0) 
through  the  first  pressure  minimum,  but  they  "separate 
out"  as  the  second  minimum  point  is  approached.  This 
clearly  indicates  that  for  an  accurate  measurement  of  the 
bulk  toss  parameter  G,  and  hence  the  admittance  of  the 
tested  solid  propellant,  the  standing  wave  structure  has  to 
be  measured  through  a  distance  covering  at  least  two 
standtog  wave  pressure  minima.  To  accomplish  this,  the 
number  of  pressure  transducers  has  been  increased  from 
ten  to  fifteen.  A  typical  distribution  of  the  pressure 
transducers  along  the  tube  is  shown  in  Fig.  4.  It  is  worth 
noting  that  the  more  accwate  piezoelectric  pressure 
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transducers  arc  located  near  pressure  minima  where  the 
greatest  experimental  accuracy  is  required. 
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Typical  admittance  and  gas  phaae  loss  data  obtained 
during  tests  conducted  with  the  an  A-IS  composite  solid 
propeUant  at  300  psig  chamber  pretaure  are  presented  in 
Pigs.  6  through  t  where  the  abscissa  “Recortf, 
corresponds  to  time.  Figure  4  shows  the  time  variation  of 
the  real  and  imaginary  parts  of  the  admittance  Y  and  the 
bulk  gas  phase  loss  coefficient  G  during  a  test  conducted 
with  an  A-lt  propellant  at  the  frequency  of  479  He.  It  is 
noted  that  the  real  part  of  the  admittance  Yg  varies  little 
with  time  while  energy  losses  in  the  gas  phase  very 
significantly  with  time.  These  figures  demonstrate  the 
ability  of  the  impedance  tube  technique  to  simultaneously 
measure  the  time  variation  of  acoustic  energy  gains  and 
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Typical  Results 

Considering  the  complexity  of  the  propellant 
admittance  measurement  technique,  the  repeatability  of 
the  experimental  data  was  of  much  concern.  The 
repeatability  achieved  is  demonstrated  in  Fig.  3  showing 
the  time  evolution  of  the  pressure  amplitudes  measured  at 
distances  of  2.0  and  11.73  inches  from  the  propellant 
surface  in  two  different  tests  conducted  with  an  A- 1 3 
propellant.  Examination  of  this  figure  indicates  remarkable 
similarity  in  the  time  evolution  of  die  amplitudes  and 
excellent  agreement  between  the  computed  values  of  Y 


o  wtanci.mm. 


ns.  s  BSftar  auns  a-issrs  miss  mo 


l 

I 


M 

ai 

-6J 

ai 

»  t»  »  s  r  x 
■tcoao 


RECORD 


nc.  A.  rrnCAL  DATA  MSASUMB  0UBMC  A  6IW  RUN 

a-i*  ttwaORM,  •  •  *r>  m,  r  •  **•  m*. 


losses  inside  the  tube  during  a  test.  Figures  7  and  S 
describe  the  frequency  dependence  of  the  real  part  of  the 
admittance  Yg  and  the  bulk  gas  phase  loss  parameter  G 
obtained  from  experiments  conducted  with  an  A* It 
propeUant.  The  curve  indicate  that  this  quantity  peaks 
at  around  1000  Hz. 

Figures  9  through  It  provide  comparisons  between 
the  predicted  (based  upon  the  determined  values  of  the 
propellant  admittance  and  the  parameters  G  and  C) 
impedance  tube  wave  structure  (i.e.,  amplitude  and  phase) 
and  the  measured  acoustic  pressure  data  at  three  different 
instances  during  a  test  conducted  with  an  A-lt  propellant 
at  the  frequency  of  023  Hs  and  chamber  pressure  of  300 
psig. 


These  figures  demonstrate  die  dependence  of  the 
Impedance  wave  structure  upon  the  gas  phase  losses,  which 
are  described  by  the  parameter  G  which  increases  in 
magnitude  from  9  to  19  during  the  test  period  covered  by 
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thaw  figures.  The  effect  of  the  change  in  G  becomes 
particularly  apparent  whan  one  examinee  Pigs.  12  and  1* 
that  dv»w  the  manner  in  which  the  slope  of  the  phase-space 
curve  changes  near  the  second  minima  point  as  the  value  of 
G  increased  from  10  to  IS  between  two  instants  of  time. 
Yet,  while  the  values  of  G  change  with  time  the  calculated 
valuee  for  the  reel  pert  of  the  admittance  Yt  remain 
virtually  the  same,  Mteating  an  almoet  constant  addition 
of  acoustic  energy  to  the  gas  phase  by  the  burning 
propellant*  Finally,  it  is  important  to  note  the  excellent 
agreement  shown  in  Pigs.  9  through  it  between  the 
calculated  wave  structure  (based  upon  the  determined 
values  of  the  propellant  admittance,  G  and  C)  and  the 
measured  expert  mimtalj  agreement  that  provide  further 
support  to  the  applicability  of  the  developed  experimental 
technique. 
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Conclusion* 

Significant  improvements  haw  bean  achieved  in  the 
data  acquisition  and  data  reduction  procedure*  in  the  solid 
propellant  impedance  tube  experiment  as  demonstrated  by 
the  results  provided  in  this  paper.  It  has  been  shown  that 
the  real  part  of  the  burning  propellant  admittance  YK 
remains  constant  during  the  steady  state  bum  period  while 
ps  phase  losses  may  very  significantly  during  the  same 
time  period.  Reproducibility  of  the  measured  quantities  has 
been  established  as  welt  as  the  capability  of  this 
experimental  technique  to  measure  gas  phase  losses  in  the 
impedance  tube  simultaneously  with  the  measurement  of 
the  burning  solid  propellant  surface  admittance.  Finally,  it 
has  been  mown  that  the  determined  propellant  admittance 
Mid  gas  phase  losses  result  In  impedance  tube  wave 
structures  that  are  in  excellent  agreement  with  the 
measured  acoustic  pressure  data. 
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Abstract  elevated  downstream  presstres  to  the  projectile  base  has 


Experimental  studies  of  base  pressure  manipulation 
for  an  axisymmetric  model  at  Mach  3  with  simulated  and 
actual  extemal/baae  burning  are  described.  Early  work 
using  contoured  test  section  wails  and  cold  gas  base  region 
injection  is  reviewed  to  demonstrate  wake  detail  and  length 
scale  changes  under  the  influence  of  simulated 
extemal/base  burning.  Tests  with  actual  combustion  of 
radmlly  and  axially  injected  hydrogen  are  then  reported. 
Outs  taming  performance  values  with  significant  base  drag 
reliction  is  diown  for  injection  end  burning  (Erectly  in  the 
near. wake  (best  burning).  Current  attempts  at  radial 
injection  and  burning  in  the  free  stream  (external  burning) 
have  not  yet  succeeded.  These  tests,  however,  have  defined 
an  envelope  within  which  external  burning,  if  feasible,  will 
presently  be  achieved. 

Nomenclature 


*b 

3 

model  base  area 

I 

s 

injection  flow  parameter  (Figure  8) 

‘m> 

X 

specific  impulse  (Figure  8) 

M 

X 

Mach  number 

fit 

s 

molecular  weight 

^Bleed 

X 

injection  mass  flow 

P 

3 

static  pressure 

•»! 

X 

free  stream  static  pressure 

Po 

3 

total  pressure  (Pitot) 

Pol 

X 

free  stream  total  pressure 

Pb 

= 

base  pressure 

Pbo 

3 

undisturbed  base  pressure 

R 

- 

model  radius 

RSP 

= 

rear  stagnation  point 

To 

X 

stagnation  temperature 

T0l 

X 

free  stream  stagnation  temperature 

V1 

X 

free  stream  velocity 

’I 

s 

free  stream  density 

Y 

X 

ratio  of  specific  heats 

Introduction 

The  application  of  burning,  either  in  the  near  wake 
(base)  or  in  the  free  stream  adjacent  to  the  near  wake 
(external),  can  provide  significant  base  drag  reduction  for 
projectiles  operating  at  air  breathing  altitudes.  The  coupling 
of  attractive  performance  values  with  hardware  simplicity 
makes  this  a  particularly  interesting  propulsion  concept. 

The  feasibility  of  external/base  burning  for 
propulsion  in  which  the  subsonic  near-wake  transmits 
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been  established.  At  present,  however,  no  analytical  model 
exists  which  adequately  treats  the  complex  base  flows  that 
result  from  combustion  in  and  near  the  wake.  A  major 
contribution  in  understand^  these  flows,  establishing  and 
optimizing  performance  criteria  and  providing  (Erection  for 
theory  can  be  supplied  experimentally  through  detailed 
pressure  and  temperature  measurements  in  the  wake  region 
for  various  modes  of  extemal/baae  combustion.  Cirrent 
studies  in  the  School  of  Aerospace  Engineering  at  Georgia 
Tech  seek  to  provide  this  information. 

A  versatile  experimental  facility  has  been 
established  for  both  simulation  and  actual  combustion 
testing  of  axisymmetric  extemal/base  burning 
configurations  at  Mach  3  with  a  5.72  cm  dia.  model  in  a 
15.2  cm  dia.  test  section  (Reynolds  number  based  on  model 
diameter  *  3.0  x  10*).  The  blow-down  facility  is  shown  in 
Figures  1  and  2.  The  hollow  cylindrical  model  is  momted  in 


Figure  1.  Extemal/Base  Burning  Test  Facility 
Schematic. 


Figure  2.  Extemal/Base  Burning  Test  Section 
Schematic. 
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the  subsonic  (low.  puses  through  the  nozzle  throat  and 
terminates  in  the  plexiglass  test  section.  Gases  (or  base 
region  injection  are  brought  in  through  (our  at  the  eight 
support  strut  access  ports.  The  remaining  tour  ports  are 
used  to  bring  model  surface  static  pressure  tubes  to  the 
external  data  acquisition  system.  Detailed  pressure  and 
temperature  wake  measurements  are  made  with  computer 


1.  D  and  HI).  One  additional  test  section  (IV)  was  (abricated 
to  show  the  eiiects  ot  discrete  burning  lor  comparison  with 
the  equivalent  axisymmetric  case  (11).  Each  test  section 
was  translated  relative  to  the  model  base  to  investigate  the 


controlled  traverses  ot  probes  introduced  through  the 
constant-area  <g((user  downstream  oi  the  test  section.  Run 
times  range  from  4-7  minutes  depending  on  the  stagnation 
pressure  selected.  Tunnel  (low  is  not  heated  and  stagnation 
temperatures  <ki(t  downward  from  ambient  to 
approximately*20°C  during  a  typical  tut. 

Protram  Review 

The  initial  phase  o(  the  present  comprehensive 
experimental  study  involved  simulation  of  various 
external/ base  burning  modes.  These  results  have  been 
reported  previously  and  will  be  summarized  here  briedy  to 
emphasize,  important  (indings  and  to  trace  program 
evolution.1’ 

First,  measurements  were  made  to  provide 
undisturbed  Mach  3  axisymmetric  base  (low  detail  (or 
theoretical  model  comparisons.  Typical  results  are  shown  in 
Figures  3  and  *  in  which  radial  and  axial  Mach  number  and 
pressure  profiles  are  presented.  These  measurements  reveal 
presence  of  a  base  recirculation  region,  wake  thickness 
variation,  streamline  curvature,  base  pressure  and  near¬ 
wake  length  as  defined  by  the  location  of  M  c  1  on  the 
centerline. 


Figure  3.  Near  Wake  Radial  Mach  Number  and  Static 
Pressure  Distributions  (Zero  External 
Compression). 


Second,  tests  with  simulated  axisymmetric  external 
burning  were  completed  using  contoured  test  section  walls. 
The  wall  contours  were  designed  to  focus  compression 
fields  on  the  stake  representing  a  reasonable  fixed  total 
heat  addition  applied  in  the  free  stream  over  three 
successively  shorter  heating  lengths  (Compression  Sections 


Figure  4.  Near  Wake  Centerline  Pressure  and  Mach 
Number  Distributions  (Zero  External 
Compression). 


importance  of  disturbance  location.  Figure  5  shows  base 
pressure  elevation  to  free  stream  values  and  above  (or 
most  cases  with  an  optimum  axial  position  depending  on  the 
simulated  heating  severity.  Base  pressure  rise  is  reduced  if 
disturbances  are  applied  discretely  rather  than 
axisymmetricaliy.  Nevertheless,  for  all  cases  examined, 
base  pressure  was  still  boosted  sufficiently  to  eliminate 
base  drag.  The  effect  of  axisymmetric  external 
compression  on  neaf-wake  length  scales  is  shown  in  Figure 
(.  Wake  lengths  are  strongly  controlled  by  external 
disturbance  location  as  reflected  by  changes  in  the  position 
of  the  rear  stagnation  and  M  =  1  points  on  the  centerline.  In 
addition,  wake  lengths  decrease  with  increasing 
compression  severity  for  a  given  axial  location. 

Further  simulation  studies  of  discrete  external 
burning  were  completed  using  solid  pegs  to  model  the 
blockage  effects  of  radial  fuel  injection  jets  (no  external 
compression).  Figure  7  compares  wake  centerline  structure 
for  this  case  with  the  original  dean  base  data.  The 
disturbances  generated  by  the  pegs  result  in  a  mild 
decrease  in  base  pressure  and  a  significant  shortening  of 
the  near-wake. 

The  third  segment  of  testing  Involved  simulation  of 
base  burning  with  cold  gas  injection  into  the  near-wake 
region.  The  model  was  fitted  with  a  porous  sintered-metal 
base  plate  for  uniform  Injection  and  tests  were  run  with 
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Figure  5.  Base  Pressure  Manipulation  with 
Axisymmetric  and  Discrete  External 
Compression. 


-three  gases  displaying  a  range  oi  molecular  weights.  Base 
pressure  elevation  and  performance  values  for  nitrogen 
(air),  helium  and  hydrogen  are  presented  in  Figure  8 
showing  the  expected  gains  with  decreasing  molecular 
weight  (increased  volumetric  flow).  Figure  9  gives  a 
successful  correlation  of  base  pressure  rise  with  injection 
rate  parameter  modified  by  a  factor  involving  molecular 
weight  and  specific  heat  ratio  of  the  injected  gas  .  Figure 
10  gives  typical  static  pressure  and  Mach  number 
distributions  in  the  near  wake  for  a  moderate  air  Injection 
rate.  A  recirculation  region  persists  lor  this  case,  although 
the  region  has  been  axially  displaced  with  base  bleed.  Tests 
at  higher  injection  rates  revealed  that  the  recirculation 
bubble  can  be  completely  blown  away. 

The  final  phase  of  non-combustion  work  combined 
simulated  axisymmetric  external  compression  with  cold  gas 
base  bleed.  Wake  length  and  base  pressure  alterations  for  a 
range  of  base  injection  parameter  values  with  and  without 
external  compression  are  shown  in  Figure  il.  Without 
compression,  base  pressures  rise  modestly  and  wake  lengths 
increase  significantly  with  Increasing  base  bleed.  When 


Figure  6.  Effects  of  External  Compression  Strength  and 
Location  on  Near  Wake  Length  Scales. 


Figure  7.  Comparison  of  Centerline  Pitot  Pressure, 
Static  Pressure  and  Mach  Number 
Distributions  for  the  Peg  and  Clean  Base. 
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Current  Work 

After  completion  of  external/base  burning 
simulation  studies,  the  facility  was  modified  for  controlled 
base  region  injection  and  burning  of  hydrogen.  Initial  tests 
were  designed  to  evaluate  performance  and  detail  wake 
structure  for  hydrogen  base  burning  with  injection  directly 
into  the  near- wake  through  a  porous  base  plate.  The  tunnel 
was  operated  at  930  kpa  (135  psia)  simulating  Mach  3  flight 
at  approximately  10,700  m  (35,000  ft).  However,  since 
tunnel  heating  is  not  available,  the  resulting  stagnation 
temperature  is  approximately  330°C  (600  F)  below  the 
corresponding  Mach  3,  10,700  m  value.  Under  these  severe 
conditions,  attempts  at  base  region  ignition  with  a  6000 
volt  A.C.  spark  source  were  unsuccessful.  As  an 

alternative  to  more  powerful  spark  sources,  a  pyrotechnic 

compound  paste  was  employed  to  achieve  ignition.  The 
paste  was  coated  on  the  surface  of  an  "L"  shaped 

cylindrical  igniter  probe  introduced  to  the  tunnel 

downstream  of  the  test  section.  The  active  portion  of  the 
igniter  was  located  along  the  wake  centerline  from  two  to 
five  model  radii  downstream  of  the  base.  Paste  firing  was 
initiated  with  an  embedded  hot  nichrome  filament.  A 


Figure  12.  Pyrotechnic  Igniter  Firing  in  Near- Wake  with 
Hydrogen  Base  Bleed. 


typical  firing  is  shown  in  Figure  12. 

Early  tests  with  burning  revealed  instability  and 
flameout  upon  removal  of  the  igniter  probe.  The  problem 
was  overcome  by  installing  four  tabs  of  catalytic  platinum 
gauze  on  the  base  plate.  The  final  igniter  configuration 
utilizing  a  consummable  base  mounted  sting  is  shown  in  the 
top  panel  of  Figure  1 3. 

Typical  operation  with  hydrogen  base  burning  at 
Mach  3  is  shown  in  Figure  14.  For  injection  rates 
greater  than  I »  4.0  x  10*’  the  visible  combustion 
region  extends  four  to  five  model  diameters 
downstream  of  the  base  plane  while  remaining 
attached  to  the  base.  For  lower  injection  flows  the 
luminous  region  diminishes  rapidly  in  size  and  intensity  and 
is  closely  confined  to  the  base.  Observation  of  igniter  flare 
particle  paths  and  pitot-static  pressure  measurements  with 
burning  indicate  that,  even  at  high  injection  rates,  a 
recirculation  region  persists  in  the  wake.  Temperature 
measurements  in  the  wake  with  platinum-platinum,  10% 
rhodium  thermocouples  show  that  combustion  exists 
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throughout  the  visible  conical  region.  These  measurements 
show  further  that  burning  extends  into  the  shear  layer. 
Temperatures  approaching  the  thermocouple  useful  limit  of 
1600°C  were  only  encountered  within  two  or  three  model 
radii  of  the  base  for  even  the  highest  injection  rates 
attempted.  This  value  is  significantly  below  the 
stoichiometric  flame  temperature  for  hydrogen  and 
indicates  the  lean  nature  of  the  wake. 

The  performance  of  hydrogen  base  burning  at  Mach 
3  is  shown  in  Figure  15  for  the  range  of  flows  and 
configurations  investigated.  The  lower  portion  of  Figure  15 
shows  that  base  pressure  rises  with  increased  hydrogen 
injection  and  approaches  free  stream  static  pressure  (near¬ 
elimination  of  base  drag).  The  diminishing  returns  of 
increased  injection,  however,  are  reflected  by  the  rapid 
decline  of  specific  impulse  values  shown  in  the  upper 
portion  of  the  figure.  Nevertheless,  a  significant  reduction 
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Figure  13.  External/Base  Burning  Injection 
Configurations. 


Figure  14.  Hydrogen  Base  Burning  at  Mach  3. 
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in  bate  drag  at  impressive  values  of  I,p  have  been 
demonstrated.  For  example,  these  tests  show  an  80% 
reduction  in  base  drag  with  Isp  -  3000  sec.  and  a  30% 
reduction  with  l*p  *  12000  sec.  For  comparison.  Figure  1ft 
presents  Mach  3  performance  for  a  hydrogen  burning 
ramjet  with  various  inlet  designs.  It  is  clear  that  base 
burning  is  strongly  competitive  with  ramjet  operation  for 
this  application. 


Further  base  pressure  studies  were  completed  by 
combining  base  burning  with  two  cases  of  simulated 
axisymmetric  external  burning  and  one  case  of  simulated 
discrete  external  burning  (use  of  contoured  test  sections 
from  original  experiments).  The  first  two  panels  (left  to 
right)  in  Figure  17  show  that  the  combined  effect  of 
simultaneous  base  and  simulated  external  burning  is  very 
nearly  the  sum  of  the  individual  contributions.  The  third 
panel  reinforces  this  finding  for  an  alternate  test  section. 
In  all  three  cases  the  base  pressure  was  boosted  above  free 
stream  static  pressure  yielding  a  condition  of  base  thrust. 
Panel  three,  in  fact,  shows  an  instance  in  which  the  base 
thrust  is  sufficiently  high  to  eliminate  all  drag  on  a  well 
designed  vehicle. 

The  last  panel  in  Figure  17  shows  results  using 
simulated  discrete  external  burning  (equivalent  to 
axisymmetric  Compression  Section  II).  Once  again,  the 
measured  base  pressure  with  base  burning  and  compression 
is  very  close  to  the  simple  sum  of  the  individual 


experimental  values.  Comparison  of  panels  three  and  four 
confirm  earlier  findings  that  the  application  of  external 
burning  in  discrete  plumes  diminishes  the  base  pressure 
boost  obtained  with  axisymmetric  burning;  however,  a 
condition  of  base  thrust  is  still  demonstrated  for  this  case. 

Tests  with  radial  injection  of  hydrogen  for 
attempted  "pure  "  external  burning  were  next  attempted.  A 
modified  base  configuration  using  six  equally  spaced 
injection  ports  recessed  in  an  annular  channel  was  used  for 
these  studies  (Figure  13).  The  channel  design  was  included 
to  provide  a  cavity  in  which  favorable  mixing  and  residence 
times  would  exist  for  ignition  and  sustained  burning.  As 
with  the  base  burning  tests,  a  pyrotechnic  compound  served 
as  the  source  of  ignition.  The  channel  was  lined  with  the 
compound  and  flared  with  a  hot  ni chrome  filament. 
Platinum  gauze  was  used  as  a  catalyst  in  the  channel  but 
not  in  the  base  region. 

Initial  tests  were  performed  with  subsonic  injection 
through  six,  9.3  mm  orifices  yielding  Jet  velocities  ranging 
from  50-170  m/sec  (I  =  3-11  x  10*4  ).  Ignition  was  readily 
achieved  for  I  >  4  x  10~*  with  burning  apparent  both 
confined  within  the  channel  and  in  the  wake  as  shown  in 
Figure  18  (wake  burning  outlined).  Reduction  of  hydrogen 
flow  resulted  in  progressively  weaker  burning  (decreasing 
luminosity)  in  the  base  region  with  wake  combustion 
terminating  abruptly  at  I*»  3.25  x  10~* .  Burning  in  the 
channel  persisted,  however,  for  even  lower  flows  and  was 
made  evident  by  the  continued  glow  of  the  platinum  gauze. 
Base  burning  probably  could  have  been  extended  to  lower 
injection  rates  by  using  platinum  mounted  on  the  base,  but 
stKh  efforts  were  not  pursued. 

Performance  values  for  "external"  burning  with 
subsonic  discrete  Jet  radial  Injection  of  hydrogen  are 
included  in  Figure  15.  As  might  be  expected  from  test 
observations,  this  mode  of  operation  more  resembles  base 
burning  with  drect  injection  rather  than  external  burning 
with  radial  injection  and  discrete  plume  combustion  in  the 
free  stream.  The  loss  of  wake  burning  is  shown  clearly  in 
the  precipitous  drop  in  base  pressure  and  corresponding 
S.p.  With  burning  in  the  channel  alone,  the  base  pressure 
is  Very  nearly  the  same  as  for  the  radial  injection  model 
configuration  with  zero  hydrogen  flow.  Since  burning  in  the 
channel  provides  a  negligible  improvement  in  base  pressure. 
Its  primary  value  must  be  considered  as  that  of  a  possible 
flamehoider  for  the  radially  injected  stream  at  some 
sacrifice  in  I(p  .  Base  pressure  and  I ,»  values  for  both  the 
base  and  subsonic  radial  injection  mooes  converge  at  the 
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Figure  17.  Mach  3  Hydrogen  Base  Burning  withSimuiated  External  Burning. 


Figure  IS.  External/ Ba»e  Burning  of  Hydrogen  et  Mach  3 

with  Subsonic  Radial  Injection. 


higher  injection  rates.  This  convergence  suggests  that,  with 
larger  flows,  increasing  percentages  of  the  radially  injected 
hydrogen  escape  the  channel  and  are  entrained  and  burned 
in  the  near-wake  region. 

The  results  of  the  subsonic  radial  injection  tests 
indicate  that,  in  order  to  have  some  chance  for  true 
external  burning,  significant  free  stream  penetration  with 
the  fuel  jets  must  be  achieved.  To  that  end,  the  original 
orifices  were  plugged  and  re-drilled  to  1mm  dia.  for  choked 
operation  at  flows  comparable  to  those  used  in  the  previous 
subsonic  Injection  tests.  Platinum  gauze  tabs  were  fastened 
to  the  base  as  well  as  in  the  channel.  Attempts  to  ignite  at 
I  *  I  x  10-4  «nd  1 1  x  10*4(jet  Mach  number^  based  on  cavity 
to  free-stream  pressure  ratio,  of  2.1*  and  2.42 
respectively)  were  repeatedly  unsuccessful.  Local  jet 
burning  was  observed  within  the  groove  during  the 
pyrotechnic  flare?  however,  once  the  source  of  ignition  was 
removed,  all  burning  ceased.  Unlike  earlier  tests  with 
subsonic  injection,  no  combustion  was  initiated  in  the 
wake  region.  It  it  apparent  that  the  higher  injection 
velocity  used  to  obtain  free  stream  penetration  for 
potential  external  burning  resulted  in  channel  conditions 
unfavorable  lor  sustaining  an  ignition  source.  In  addition, 
the  absence  of  bate  combustion  Indicates  either  that  the 
high  velocity  injection  has  altered  the  mechanism  by  which 
channel  ighition  begets  base  ignition  or  that  the  increased 
Ifrss-strsam  jet  penetration  has  significantly  reduced 
hydrogen  entrainment  in  the  stake.  The  need  for  testing 
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radial  infection  velocities  between  the  extremes  Just 
described  is  evident.  These  tests,  in  which  orifice 
diameters  are  enlarged  incrementally  to  provide  the 
desired  velocity  range  addle  maintaining  comparable 
hydrogen  flows,  are  in  progress.  Included  in  these  tests  are 
attempts  at  igniter  firings  in  the  base  region  as  well  as  in 
the  channel  to  ascertain  the  presence  of  hydrogen  in  the 
wake. 

Concurrent  with  tests  varying  radial  injection  Mach 
number  are  those  in  which  supplemental  base  bleed  is 
introduced.  The  original  radial  injection  base  has  been 
modified  for  these  experiments  as  shown  in  Figure  13.  One 
test  series  with  this  configuration  using  only  channel 
ignition  has  been  completed.  Hydrogen  injection  was  spilt 
approximately  60%  -  S0%  in  favor  of  base  injection  with 
attempted  ignition  at  I  =  10.9  x  I0"4  (based  on  total 
hydrogen  flow).  This  condition  resulted  in  a  radial  Jet  Mach 
number  of  2.07  and  sufficient  hydrogen  base  bleed,  as 
determined  from  original  base  burning  tests,  to  assure  wake 
combustion.  Sustained  channel  burning  was  achieved  for 
this  case  but  base  burning  was  not  initiated.  Totalhydrogen  2. 
injection  was  decreased  in  steps  to  1  a  3.67  x  10*4  (M  a 
1.30)  with  no  loss  in  channel  burning.  Base  pressure 
remained  nearly  constant  throughout  the  range  and  was 
only  about  20%  higher  than  for  this  identical  base 
configuration  with  zero  injection.  These  results  confirm 
the  earlier  suggestions  that  the  higher  radial  injection  rates  3. 
tend  to  destroy  the  mechanism  which  induces  wake  ignition 
from  channel  ignition  and  that  channel  burning  has  little 
effect  on  base  pressure.  It  is  probable  that  most  of  the 
modest  base  pressure  rise  shown  is  attributable  to  the 
hydrogen  bleed  directly  into  the  base  region.  The  most 
significant  result  of  this  test  series,  however,  is  the  lack  of 
sustained  burning  in  the  radial  Jett  away  from  the  channel 
confines  even  though  a  continuous,  stable  soured  of  ignition 
is  present. 

The  experiments  performed  to  this  time  have 
defined  the  limits  within  which  external  burning  with 
hydrogen  can  be  obtained  with  the  present  facility.  At  one 
extreme,  low  subsonic  radial  injection  results  in  an 
essentially  base  burning  condition.  At  the  other  extreme, 
injection  at  Mach  numbers  exceeding  about  2.1  renders 
maintenance  of  a  flame  holding  region  impossible  for  this 
configuration.  If  true  external  burning  is  to  be  achieved, 
it  will  exist  within  the  envelope  of  injection  Mach  numbers 
and  mass  flows  under  current  examination. 

Summary  and  Conclusions 

Tests  with  axial  injection  and  burning  of  hydrogen  in 
the  near-wake  have  shown  that  significant  base  drag 
reduction  and  excellent  performance  can  be  combined. 
Drag  reductions  of  80%  and  30%  have  been  demonstrated 
tor  specific  impulse  values  of  3000  sec  and  12000  sec 
respectively.  Near-elimination  of  base  drag  appears  to  be 
an  upper  limit  for  this  mode  of  external  bunting,  however, 
and  can  only  be  achieved  with  heavy  penalties  in  Icp  . 

Pure  external  burning  with  radial  injection  and 
combustion  in  the  free  stream  offers  the  hope,  of 
generating  base  thrust  by  elevating  wake  pressures  above 
tree  stream  levels.  Tests  with  low  radial  injection 
velocities  (M  <  0.1)  result  in  no  appreciable  free  stream  Jet 
penetration.  Hydrogen  is  apparently  swept  into  the  wake 
.  and  bums  as  though  it  has  been  injected  directly  through 
the  base.  Accordingly,  this  mode  performs  similarly  to  but 
less  efficiently  than  the  drect  base  burning  mode.  Radtol 
injection  with  comparable  mass  flows  through  choked 
nozzles  to  increase Jet  penetration  has  not  yet  yielded  true 
external  burning.  Flameholdar  combustion  for  jet  Mach 
numbers  greater  than  about  2.1  cannot  be  sustained. 
Successful  flameholder  operation  has  been  demonstrated 


for  jet  Mach  numbers  from  1.3  to  2.1,  but  propagation  of 
combustion  into  the  free  stream  is  not  evident.  Operation 
in  this  mode  gives  base  pressures  only  slightly  higher  from 
those  with  no  injection  at  all. 

Current  tests  exploring  jet  injection  Mach  numbers 
between  the  extremes  so  far  investigated  should  reveal 
whether  or  not  pure  external  burning  is  feasible.  These 
tests  are  being  run  with  base  as  well  as  flameholder 
ignition  to  ascertain  the  presence  of  combustible  mixtures 
in  the  wake.  In  addition,  cases  are  being  examined  in  which 
direct  base  bleed  is  combined  with  radial  injection  to  verify 
and  supplement  the  results  already  in  hand. 
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Abstract 

Theoretical  acoustics  applied  to  the  pipe  flow  turbulence, 
problem  in  a  configuration  which  was  experimentally  investigated  yields . 
an  analytical  formula  for  random  pressure  fluctuations  generated  jby  the ' 
turbulence.  In  order  to  calculate  the  pressure  fluctuations,  it  is  found  j 
that  certain  spatial  correlations  of  the  fluctuation  in  the  square  of  the  i 
axial  velocity  are  needed.  These  are  measured  by  hot  film  anemoinetry,  j 
as  are  the  pressure  fluctuations,  in  a  fully  developed  pipe  flow  and  the j 
theory  is  shown  to  be  semi-quantitative;  a  lower  bound  to  the  pressure  j 
fluctuation  emerges.  Application  to  the  problem  of  pressure  fluctuation 
in  rocket  motors  and  attendant  vibration  are  indicated.  1 

■  •  •  •  .  if»t  •  S  •  *  .t  u  s  -  -  .*1 

Nomenclature  j 

=  pipe  cross  section  area  i 

=  pipe  radius  i 

a  speed  of  sound 

a  quantity  defined  by  Eq.  (9)  j 

a  Green's  function  for  H  j 

a  Green's  function  for  plane  wave  acoustics  problem 

a  Bernoulli  enthalpy  1 

a  Bessel  function  of  order  m 

a  duct  length 

a  outward  unit  normal  vector  i 

a  pressure 

a  cross  spectral  density  between  signal  i  and  signal  j  j 

a  Eigenvalue  in  the  problem  for  H  • 

a  entropy 

a  time 

a  sample  time  in  Fourier  transform 

a  vortical  velocity  vector 

a  volume  I 

a  velocity  vector 

a  potential  velocity  vector..  1  1 

a  Cartesian  coordinate  in  ixn  direction  ! 

a  distance  along  duct  1 

a  distance  along  duct  I 

j 

a  quantity  defined  by  Eq.  (10) 

a  Dirac  delta  function  1 

*  density 

a  radian  frequency 

a  velocity  potential 
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iSubscripts 

transverse  part 
plane  approximation 

correlation  scale  I 

value  at  infinity  in  the  potential  field 
Fourier  transform  i 

x- direction  j 

Supercripts 

'  =  fluctuation 

-  =  mean  or  ensemble  average 

|  Introduction  | 

Pressure  fluctuations  inside  of  a  rocket  motor  cavity  abe  one 
source  of  motor  vibrations.  Typically,  the  broad  band  random  pressure 
fluctuations  which  are  normally  considered  reasonable  to  tolerate  are  of 
the  order  of  j_%_rins  magyjMdetWhen^par,^  with  the  rneanchambe£ 
pressure.  Depending  upon  the  frequency  content  of  this  background 
noise  in  the  pressure,  however,  this  1%  level  can  cause  vibration 
troubles.  Especially  if  die  excitation  frequencies  overlap  a  longitudinal 
mode  frequency  of  the  chamber  cavity,  strong  longitudinal  resonances 
may  occur  with  attendant  thrust  fluctuations.  Such  fluctuations  in 
thrust  can  be  more  than  the  1%  level  of  the  pressure  fluctuations 
because  the  pressure  field  is  spatially  variant  and  differentially  ^cts  on 
the  head  and  aft  ends  of  the  cavity.. This  kind  of  problem  is  currently 
believed  to  exist  in  one  large  system.  ' 

There  are  several  potential  sources  for  the  background  pressure 
fluctuations,  and  many  of  these  have  been  considered  in  R4f.  (2). 
However,  the  results  of  Ref.  (2)  contain  an  ambiguity  in  interpretation 
due  to  the  simplified  theory  employed,  and  it  is  desireable  to  use  a  more 
accurate  theory  to  treat  the  problem.  The  approach  here  is  to  look  at 
the  problem  of  turbulence  induced  pressure  fluctuations,  both 
theoretically  and  experimentally,  with  a  view  towards  prediction'  of  the 
pressure  fluctuation  level  in  rocket  motors,  due  to  the  turbulenc^  field. 
There  will  necessarily  be  left  out,  at  Xtus  point,  other  potential  sources 
of  vibration  such  as  combustion  noise.1  '  I 

The  approach  here  is  to  first  look  at  a  relatively  well 
bharacterized  turbulence  field,  that  of  pipe  flow,  encountering  a 
focket-like  flow  termination,  that  of  a  choked  de Laval  nozzlfe.  The 
purpose  is  to  see  the  accuracy  with  yhich  the  pressure  fluctuation 
levels  may  be  predicted  by  theory,  given  appropriate  measurements  of 
the  turbulence  field.  If  satisfactory  answers  may  be  obtained,  the 
theory  may  be  applied  to  rocket  motors  as  further  information  becomes 
available  on  the  rocket  motor  turbulence  field,  as  given,  for  example,  j 


t 

oo 

tor 


to 

1 


by  the  computations  of  Ref.  (3).  For  the  moment,  however,  this  may  be, 
regarded  as  a  fundamental  investigation  of  turbulence  generated^ 
pressure  fluctuations  in  a  given  type  of  interior  flow  problem. 


Analysis 


Consider  the  configuration  of  Fig.  1,  where  a  pipe  flow  is 
terminated  by  a  choked  nozzle  and  is  also  the  configuration  in  the 
experimental  results  to  follow.  Except  very  near  the  walls  the  inlet  flow 
is  nearly  potential;  there  follows  a  transition  to  turbulent  flow  and 
finally  a  fully  developed  purely  vortical  turbulent  pipe  flow.  In  the  fully 
developed  turbulent  region  it  is  the  mean  flow  which  is  purely  vortical,, 
but  superimposed  is  a  turbulent  fluctuation  which  has  pressure 
fluctuations  associated  with  it.  The  calculation  of  the  pressure  field  isj 
the  issue  here. 


Consider  the  formulation  of  the  conservation  equations  as  given ; 
by  Yates  and  Sandri. 
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(5) 


In  the  above,  all  effects  of  transport  phenomena  are  neglected  so  that 
when  considering  fluctuations  it  is  understood  that  the  frequency, 
domain  under  consideration  is  that  of  the  energy  containing  eddies  of, 
the  turbulent  flow.  A  companion  equation  to  Eq.  (2)  has  been  omitted;  it! 
is  the  vorticity  conservation  equatibn.  H  is  the  Bernoulli  enthalpy  given ; 
by 


H  -  h  ♦ 


&  +  IzrJ — 
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and  plays  a  fundamental  role  in  the  generation  of  noise  in  Eq.  (1).  In  thej 
potential  flow  region  of  the  pipe  flow  "2H  =  0  and  H  =  H  a  constant.j 
In  the  fully  developed  pipe  flow  region  such  a  relation  does  not  hold,  but 
Eq.  (2),  using  the  consequence  of  Eq.  (9),  becomes 
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It  will  be  assumed  that  the  majority  of  the  experimental  configuration 
is  filled  with  the  fully  developed  pipe  flow  so  that  Eq.  (7)  is  to  be 
solved. 

It  is  assumed  that  the  potential  velocities  are  much  less  than  the 
vortical  velocities,  which  may  be  verified  after  the  fact,  so  that  Eq.  (7)f 
becomes 
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The  genesis  of  Eq.  (2)  is  the  momentum  equation 

VH  -  V  (  -  V  .  7  V 

which  is  zero  on  the  walls  of  the  tube  because  all  velocities  must 
vanish.  At  the  head  end  of  the  tube  7H  =  0  since  the  flow  is  a  potential 
flow.  At  the  nozzle  end,  however,  more  care  is  required.  What  is 
needed  to  solve  Eq.  (8)  is  the  boundary  condition  7  H*n  on  all 
boundaries  of  the  flow  region.  j 

! 

Splitting  the  velocity  into  its  mean  plus  fluctuating  components! 
and  noting  that  at  the  nozzle  end  ^=0  and  V  =  uj  (r)  ,  the  xj  =  x! 
component  of  Eq.  (9)  becomes 
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In  Eq.  (10)  the  fluctuations  in  the  potential  field  have  been  neglected  as; 
compared  with  those  in  the  vortical  (turbulence)  field.  j 

Now  it  is  presumed  tht  the  vortical  field  is  known  or  can  bej 
measured  so  that  the  right  hand  sides  of  Eqs.  (8)  and  (10)  are  known.; 
The  solution  may  be  written  in  terms  of  the  Green's  function  for  the 
Poisson  equation 
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7g  •  n  •  0  on  the  side  wall  and  exit  plane 
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g  •  0  at  the  entrance  plane 
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whereby 
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The  GreenS  function  may  be  constructed  by  standard  methods 


and  is  given  by 
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where  satisfies 
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and  s _ is  determined  by 

mn 


3m  (s  a)  3  o 
mn  m  mn 


The  term  in  the  solution  with  m  =  n  =  0  has  only  an  axial  dependence! 
and  corresponds  to  cross  section  averages  of  quantities.  Emphasizing] 
this. 
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where  the  t  subscript  denotes  the  part  which  is  due  to 
dimer. sions  of  the  duct. 


trail: 


averse 


The  crucial  step  is  now  taken  to  ignore  the  higher  terms  in  gj 
other  than  goo  .  It  will  be  argued  later  that  this  will  yield  a  loweij  bound] 
on  the  pressure  generated  by  the  turbulence.  Constructing  gfromEqs. 
(1*),  it  is  found 
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The  fluctuation  in  H,  H*  =  H^-  H  ,  is  then  constructed  from  Eq.  <12* 
after  considerable  manipulation  and  noting  that  v.  ij  =  0.  The  result  is  I 


»'  -  i  Jd*  fry.,)' 


(15V 


This  result,  remarkable  in  its  simplicity,  says  that  Bernoulli  enthalpy 
fluctuations  are  merely  created  by  the  fluctuations  in  the  square  of  the 
axial  velocity  components.  I  ! 

i 

Now  moving  to  construction  of  the  unsteady  potential  field,  it  isj 
readily  shown  that  for  low  Mach  number,  which  is  assumed,  ana 
sufficiently  small  fluctuations  in  the  potential  field,  which  has  Already^ 
been  assumed,  Eq.  (1)  reduces  to 
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wave! 


Since  H  is  purely  planar  by  the  above  developments,  a  plane  wavi 
solution  to  Eq.  (16)  is  adequate.  After  taking  the  Fourier  transform  ofj 
Eq.  (16),  the  solution  is  r ,  J 


or 
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where  is  the  plane  wave  Green's  function  for  the  duct.  It,  again,  mays 
be  constructed  by  standard  means  and  depends  upon  the  impedance  of 
the  duct  walls  and  end  planes.  If,  as  in-  the  experiment,  it  may  be 
assumed  that  the  walls  are  hard,  the  open  end  is  likegan  open  endea 
organ  pipe  and  the  nozzle  end  behaves  like  a  hard  wall,  an  especially] 


simple  form  for  gM  emerges.  It  is 
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Since  all  terminations  and  walls  were  assumed  perfectly  reflect  l  g>  EqJ 
(18)  contains  resonances  of  infinite  magnitude.  It  is  not  hard  to  show, 
however,  that  Eq.  (18)  is  valid  around  the  troughs  even  when  the  walla 
are  nearly  perfectly  reflecting.  Equation  (18)  will  be  used,  therefore, 
away  from  the  resonance  frequencies. 
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Finally,  the  linearized,  low  Mach  number  form  of  the  equation  of , 
state,  Eq.  (5),  using  Eq.  (6)  becomes 

p'  •  »o(H'  -  fl  >  1 
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Using  this  with  Eqs.  (15)  and  (17),  the  final  relation  for  the  pressure 
transform  becomes 

i 
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« 

a  I  +  II  (19)| 
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The  pressure  fluctuation  therefore  consists  of  two  superimposed!  parts. 

1  is  a  local  noise  due  to  the  local  turbulence  field.  II  is  a  propagational ; 
sound,  because  of  the  appearance  of  gtf  ,  and  is  the  sound  pumped  up  by 
the  entire  volume  of  turbulence.  The  frequency  content  of  I  is  the 
frequency  of  the  the  square  of  the  axial  velocity  fluctuation,  and  that ! 
of  II  is  modified  by  g  . 

Now  consider  the  construction  of  the  cross  spectral  density  of  a 
pressure  at  two  different  points,  Xj  and  Xg  •  Using  Eq.  (19), 

S12  “  ^(Xj)  I(x2)  +  SI(Xl)  II(x2)  +  SII(x1)I(x2) -,‘SII(s1)II(x2) 

Using  the  notation  T  (x)  =  (u^u^  V  ,  it  is  seen 


sI(x1)  l(x2) 
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nut  if  xj  and  Xj  are  sufficiently  separated  there  is  expected  no 
correlation  between  T  (Xj)  and  T  (^  )  soSwxl)  i(X2)«/).On  the  other 
hand,  if  xj  -  x2  so  that  an  autospictrum  of  the  pressure  is  considered 
Sjj  is  nonzero  and  is  given  by 

*11-  J«lA(r,e)  JdA(7,'§)  T  *(x,r,Q)  ^(x/r/g) 


Investigation  of  the  magnitude  of  Eq.  (20)  reveals  that  it  is  expected  to! 
be  of  order  I 
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bimuar  manipulations  on  sn(xi)li(x2)  yieia  tnat  even  wnen 
is  and  Xj  are  separated  there  is  correlation  between  the  signais  due  to 
tne  fact  mat  II  corresponds  to  propagational  sound.  In  fact  ! 
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The  point  is  that  both  x0  and  x0  run  from  0  to  t  and  where  3^3?  the  T; 
values  are  correlated.  Careful  examination  yields  the  order  of 
magnitude  result 
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so  an  autospectrum  is  being  considered,  and  the  ratio 
is 
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cos  k’C  A, 


The  exact  nature  of  this  ratio  depends  upon  the  behavior  of  the; 
correlation  area  and  volume  with  frequency,  but  generally  it  is  expected 
that  a)  the  propagational  sound  may  dominate  near  resonance  peaks 
(cos  hi  os  0)  and  b)  away  from  resonance  s  will  dominate  since  hi 
is  of  order  unity  or  larger  but  k  V_or /A*  r  is  roughly  like  j 
k  l„0~  «  mt/cw  (uja/u)  (u/c)  <<  1  for  low  Mach  number  flows  with 
turbulence  centered  about  frequencies  w  a/u  «  1 .  It  may  further  be , 
shown  that  the  cross  terms  Si  xi  and  Sn  x  are  small  compared  with  i 
either  Si  x  or  Sn  n  Consequently,  the  situation  is  as  follows:  a)  If  x^ 
and  x2  are  widely  separated  only  Sit  iishould  be  sensed  by  a  cross : 
spectrum  of  two  pressure  transducers,  b)  If  =  x2  the  auto  spectrum 
of  the  pressure  will  contain  mainly  St  t  with  contamination  by  Sim  at 
resonant  frequencies  and  c)  If  enough  information  Is  available  on  the 
spatial  correlation  of  Tw  the  pressure  should  be  calculable.  . 

It  is  argued  that  the  pressure  computed  from  Eq.  (18)  is  a  lower ' 
bound  to  the  actual  pressure.  This  conjecture  is  due  to  the  retention  of 
only  one  term,  the  Plane  approximation,  in  H  *  .  Actually,  the  Green's 
function  of  Eq.  (13)  contains  an  infinite  number  of  terms.  However, 
each  term  is  somewhat  uncorrelated  with  every  other  term  becaiuse  of 
different  transverse  function  weighting  in  the  integration  operation  of 


Mm 


mMm 


Eq.  (12).  Consequently,  the  mean  square  pressure  should  consist  of  a 
sum  of  squared  terms,  the  first  of  which  is  given  by  the  solution  asj 
given  here.  Further  work  in  completion  of  such  a  solution  to  include! 
transverse  dependence  of  g  should  be  carried  out. 

i 

Experimental  '  j 

The  tube  of  Fig.  1  is  20  ft  long,  2  inches  in  diameter,  andj 
instrumented  with  flush  mounted  microphones  and  hod  film 
anemometers.  Because  (uj  u^  Y  is  desired  from  the  hot  films,  thej 
individual  signals  are  first  passed  through  a  linearizer,  then  squared  and* 
then  filtered  to  remove  the  d.c.  part.  Ail  signals  are  recorded  on| 
magnetic  tape  for  later  Fourier  .analysis.  j  ! 

The  nozzle  for  most  runs  has  a  contraction  ratio  to  give  an: 
average  Mach  number  of  0.1  at  the  entrance  plane.  To  insure  fully! 
developed  pipe  flow  by  this  station  measurements  were  made  -of  the* 
mean  velocity  and  axial  turbulence  intensity.  These  measurements  are! 
shown  in  Fig.  2  and  are  compared  with  LauferV3'  pipe  flow  data.  Thej 
axial  extent  of  the  fully  developed  flow  was  not  determined.  { 

•  •  ►; 

First  consider  the  pressure  auto  spectrum  shown  in  Fig.  3.  Thej 

pressure  is  here  being  measured  2  inches  upstream  of  the  nozzle.  Thej 
spectrum  consists  of  a  broad  band  noise  on  which  there  are1 
superimposed  weak  resonance  peaks  at  the  1/4  wave  resonance  points! 
(here  15,  40,  70;  etc.  Hz).  This  is  the  situation  expected  from  the  theory! 
above.  Figure  4  displays  the  cross  spectral  density  magnitude  be  tween  j 
two  microphones  placed  8  inches  apart.  Clearly  what  has  been  do0e  is  to 
suppress  Sj,  .  t  (xrl  *  as  expected,  and  to  reveal  only  the 

propagational  sound/  Tnese  observations  are  in  accord  with  the 
qualitative  expectations  from  theory;  the  question  now  is  the 
quantitative  agreement.  j 

To  see  if  the  pressure  is  calculable,  what  is  needed  for  Sji  is 
given  in  Eq.  (20).  Using  the  cylindrical  coordinate  system  of  Fig.  5,  this 
becomes 


2a  tp2(p,r,0) 

H"*  P  ST(t.e)I(r,0.p.») 

°  ^(p.r.Q) 
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While  axial  symmetry  helps  somewhat,  this  still  requires  an  enormous: 
number  of  cross  spectral  densities  of  space  separated  hot  films  toj 
obtain  accurate  numbers.  Accordingly,  an  approximate  procedure  is? 
used.  A  representative  r,  point  is  taken  as  r  =  a/2  0  -  n  /2  and  a  fixed? 
anemometer  is  placed  at  this  point.  A  traverse  of  separation  distance 
p  =  0  to  a/2  is  taken  with  9.  =rt  /2  and  then  a  traverse  p  =  0  to  3  a/2  is 
made  with  9  =  3n  /2.  S_  .  is  then  approximated  by  numerical, 
integration  of  the  results.  I  j 

a/2  3a/2  ! 

-  *11  ~  T  2"  J  +  J  sn\°Xn  | 

Shown  in  Fig.  6  is  the  behavior  of  Sj?  as  a  function  of  separation 
distance ,  p  ,  for  various  frequencies.  As  expected,  the  size  scale  over  j 
which  neighboring  anemometers  are  correlated  decreases  w|th  an! 
increase  in  frequency.  Sn  is  then  calculated  by  integration  ahd  thej 
result  is  shown  on  Fig.  3.  It  is  seen  that  a)  indeed  a  lower  bound  isj 
produced  for  the  broadband  noise,  b)  the  frequency  content  is  roughly: 
correct  and  c)  this  lower  bound  is  about;  10  dB  lower  than  the, 
experimental  curve  (that  is,  the  theory  under  estimates  the  actual 
linear  rms  pressure  by  about  a  factor  of  3). 


Sn  xi  is  more  difficult  to  come  by  because  there  is  an  axial, 
integration  also  involved  in  Eq.  (21).  Looking  at  the  auto  spectrum,: 
letting  Xq  =  Xo  ♦  r  ,  and  assuming  the  Green's  function  varies  slowly ! ' 
over  the  distance  scale  where  the  axial  correlation  of  T  drops  to  zero,j 
there  results  the  approximation 


SUH  "  Po  fc4  Jdxo  8«<x»xo)8«#(x»xo)J:dC  JdA<r»0>  X 

JdA(p ,9)  ST(Xo>r>e)T(Xo+ 

The  difficulty  here  is  that  now  has  a  phase,  dependent  on  Q  , 

because  of  convection  of  disturbances,  rlowever  Sim  is  real  so  that  the 
integration  must  cancel  out  all  phase.  Consequently,  only  the  real  part 
of  Sipf.  is  used  above.  A  typical  example  of  St/~  -  a\rrx  +  r  r 
is  shown  in  Fig.  7.  The  above  formula  is  estimated  £$*  *  *0  *  » 


W5*4(^A) 


dx 


cor 


(22) 
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where 


cor 


T(x0,r,e)T(xo+C,  r,  0) 


dC 


and  the  integral  is  numerically  computed  from  the  data.  This  is  fclearlyj 
a  rather  coarse  estimate,  and,  recall,  the  Green's  function  is  noti 
accurate  near  the  resonance  peaks.  The  calculation  of  Eq.  (22)  id  shown; 
near  the  troughs  in  g^  on  Fig.  3.  Again,  a  lower  bound  is  evident  and 
the  troughs,  as  expected,  are  significantly  down  from  the  lower  bound 
estimate  of  Sjj  . 

In  the  experimental  configuration,  the  Sjjterm  is  dominating  the 
pressure  fluctuation.  Integration  of  the  spectrum  to  yield  the'  mean 
square  pressure  is  made  difficult  because  of  the  ACor  term  in  the  order 
of  magnitude  estimate  below  Eq.  (20)  .  If  Aeor  were  independent  of 
frequency  j 


<P 


•2 


5* 


would  result.  However,  since  Acor  may  be  expected  to  fall  withj 
frequency,  a  dependence  on  u  to  somewhat  less  than  the  fourth  power 
would  be  expected.  Checking  this  with  nozzles  of  differing  contraction  j 
ratio,  the  result  is  shown  in  Table  1.  j  j 


Table  1 

Scaling  of  Pressure  Level  with  Flow  Velocity 


Mass  Weighted 

Sound  Pressure  .  , 

Flow  Velocity 

Level  (dB  re  2x10°  n/mz) 

9 

« 

( 

(ft/sec) 

112 

108.3 

162 

111.8 

214 

113.8  j 

i 

The  result  is  that  <p'2  >  «  u3  approximately. 

*  * 

Concluding  Remarks  i 

What  has  now  been  gained  is  an  understanding  of  the  pressure  j 
fluctuation  generation  by  turbulence  in  an  interior  flow.  This 
understamfing  is  semi-quantitative  in  the  sense  that  a  lower  bound  to 
the  pressure  field  may  be  calculated  from  knowledge  of  ftpatial! 


... 


l 

1 


correlations  of  the  turbulence  velocities.  The  next  step  would  be  to  use 
information  on  rocket  motor  turbulence,  if  available,  together  with  an 
accurate  calculation  of  g  as,  for  example,  in  Ref.  (2),  and  to  use  these! 
pieces  of  information  to  calculate  pressure  fluctuations  in  a  rocket 
motor  cavity.  The  problem  is  that  the  necessary  correlations  are  not; 
available.  While  progress  is  being  made  in  calculation  of  the  kinetic, 
energy  of  turbulence197  the  spatial  correlations  necessary  here  do  not 
appear  to  be  forthcoming.  Consequently,  current  work  is  centering  on 
measurement  of  the  necessary  turbulence  quantities  in  a  tube  as  In  Fig.  j 
1  but  with  blowing  from  the  side  walls,  to  simulate  a  solid  rocket; 
cavity.  • 
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Figure  Captions 

Fig.  1.  Schematic  drawing  of  the  experimental  configuration. 

Fig.  2.  Comparison  of  Laufer*s  pipe  flow  data  with  that  obtained' 
in  the  current  experimental  configuration.  ! 

Fig.  3.  Measured  pressure  spectrum  at  the  wail  and  at  the  kozzlej 
entrance  plane.  Calculated  contributions  to  the  pressure  from  the 
velocity  fluctuation  data.  j 

| 

Fig.  4.  Cross  spectral  density  between  two  microjihones 
separated  by  eight  indies.  | 

Fig.  3.  Coordinate  system  for  spacewise  cross  correlations  of 
velocity  fluctuations  and  actual  experimental  points. 

Fig.  6.  Cross  spectral  density  of  transversely  separated  hot  film 
anemometers  for  various  frequencies. 

Fig.  7.  Real  part  of  cross  spectral  densities  of  axially  separate*, 
hot  film  anemometers  for  various  frequencies. 
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